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I INTRODUCTION 

t The  retrieval  and  Interpretation  of  information  from  remote  sensing 

of  the  real  nonhomogeneous  atmosphere  requires  a combination  of  accurate 
radiometric  measurements  and  reliable  computations  of  compatible  atmo- 

* 

spheric  transmittances . For  repetitive  applications,  it  is  not  practi- 
cable to  perform  complete,  rigorous  computations  of  actual  atmospheric 
transmittances  on  the  basis  of  the  characteristics  of  all  absorption  lines 
at  all  distributions  of  absorber  amounts,  pressures,  and  temperatures. 
However,  it  is  practicable  to  design  simplified,  comprehensive  transmit- 
tance representations  that  are  based  on  transmittances  computed  rigor- 
ously for  select  idealized  conditions.  Once  such  flexible  representations 
are  achieved,  they  provide  the  basis  both  for  the  interpretation  of  mea- 
surements and  for  an  a priori  simulation  of  multispectral  measurements. 

The  simulation  capability  is  particularly  important  for  an  experimental 
design  involving  passive  multispectral  radiometric  probing  as  a means 
for  acquiring  information. 

To  be  useful,  a numerical  model  must  agree  reasonably  well  with  ob- 
servations of  the  actual  atmosphere.  Unfortunately,  some  of  the  atmo- 
spheric targets  (e.g.,  clouds)  for  which  information  is  desired  through 
remote  sensing  impose  difficulties  in  data  interpretation  because  of  un- 
certainties in  their  particulate  distributions  and  radiative  properties. 
Therefore,  a logical  first  step  in  the  development  of  tools  for  infrared 
remote  sensing  studies  is  the  design  of  a suitable  transmittance -radiance 
model  for  the  molecular  (clear)  atmosphere.  Such  a model  enables  the  de- 
scription of  clear-sky  background  radiances  and  permits  Inferences  of  the 
effects  of  clouds  and  other  aerosols  on  radiance  measurements  from  com- 
parisons between  measurements  and  computations.  As  more  Information  be- 
comes available,  the  transmittance  model  can  be  altered  to  Include  partic- 
ulate constituents  of  the  atmosphere. 
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Practical  methods  for  defining  multispectral  transmittance  represen- 
tations for  arbitrary  absorber-pressure-temperature  configurations  were 

1 * 

developed  in  previous  studies,1  and  form  a basis  for  this  study. 


* 


References  are  listed  at  the  end  of  the  report. 
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II  OBJECTIVES 


► 


The  major  objective  of  this  study  was  the  establishment  of  the  basis 
for  examining  the  relative  responsiveness  of  multispectral  infrared  and 
microwave  satellite  sensors  that  can  be  applied  for  probing  of  that  por- 
tion of  the  atmosphere  (troposphere)  that  contains  clouds.  Multispectral 
measurements  are  often  used  to  generate  information  regarding  depth  dis- 
tribution of  temperature  or  absorbers  within  the  atmosphere.  At  distinctly 
different  spectral  regions  the  radiative  properties  of  a target  of  interest 
(e.g.,  cloud)  may  vary  significantly.  Exploitation  of  the  latter  charac- 
teristics might  best  be  accomplished  by  remote  sensing  in  widely  separated 
wavenumber  intervals  that  offer  the  same  response  to  the  cloud -free  atmo- 
sphere. Accordingly,  this  study  emphasized  examination  of  transmittance 
weighting  function  profiles  for  spectral  regions  responding  to  the  middle 
and  lower  troposphere  and  typical  of  current  satellite  sounder  instruments. 
The  range  of  spectral  regions  considered  extended  from  the  near -infrared 
(4  urn)  through  several  window  regions  and  bands  of  carbon  dioxide  and  water 
vapor  into  the  rotational  band  of  water  vapor  and  out  to  the  60 -GHz  region 
of  oxygen  absorption,  as  well  as  the  183-GHz  and  22-GHz  lines  of  water 
vapor. 

From  the  point  of  view  of  cloud  probing  potential,  responsiveness  to 
the  lower  troposphere  is  important.  The  broad  nature  of  typical  response 
(weighting)  functions  means  that  emission  from  the  lower  boundary  is  in- 
cluded in  the  response.  Therefore,  any  background  study  such  as  this  must 
include  variations  in  the  lower  boundary  as  well  as  those  of  the  atmosphere. 
Typical  atmospheric  variations  that  influence  the  response  include 
temperature -moisture  structural  changes  and  changes  in  viewing  angle.  Of 
course,  cloudiness  itself  would  change  with  view  angle. 

The  effects  of  clouds  on  the  radiometric  response  can  only  be  exam- 
ined properly  with  adequate  radiation  models  of  the  clouds  themselves. 
However,  in  addition  to  background  response  variations,  some  gross  cloud 
effects  can  be  estimated  from  data  summaries  compiled  in  this  study. 
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Ill  APPROACH 


Figure  1 Is  a summary  of  the  major  tasks  undertaken  to  meet  the  ob- 
jectives discussed  above.  Some  of  the  pertinent  details  of  each  step  will 
be  discussed  in  later  sections.  Both  microwave  and  infrared  radiation  are 
included  in  each  study  phase.  A summary  of  computer  programs  used  in  this 
study  is  presented  in  Appendix  A. 

Since  subsequent  transmittance  models  for  spectral  bands  should  be 
consistent  with  line-by-line  calculations  for  any  prescribed  response 
function,  it  is  necessary  to  base  the  modeling  on  the  individual  absorption 
line  parameters  (position,  strength,  half-width,  and  energy  of  the  lower 
state).  The  initial,  and  most  costly,  study  phase  is  the  use  of  the  ab- 
sorption line  parameters  in  conjunction  with  a suitable  line  shape  to  gen- 
erate detailed  files  of  monochromatic  absorption  coefficients  at  a variety 
of  pressures  and  temperatures  characteristic  of  the  atmosphere.  All  per- 
tinent molecular  absorbers  in  the  atmosphere  must  be  included  in  the  com- 
putations. The  maximum  spectral  spacing  for  the  computations  is  governed 
by  the  line  shape  and  the  pressure-dependent  half-width  at  half-maximum. 

In  addition  to  the  transmittance  model  input,  computed  absorption  coeffi- 
cients can  assist  in  the  selection  of  spectral  intervals  of  particular 
interest  or  in  the  evaluation  of  relative  performance  of  numerical  trans- 
mittance models. 

The  next  major  step  is  the  infrared  transmittance  modeling  itself. 

This  begins  with  the  calculation  of  monochromatic  transmittances  from  the 
calculated  absorption  coefficients.  A tremendous  data  compression  is  then 
obtained  by  averaging  the  transmittances  over  specified  spectral  intervals 
of  interest.  During  the  averaging,  any  prescribed  spectral  response  func- 
tion can  be  convolved  with  the  monochromatic  transmittances.  Next,  a 
selected  one-parameter  transmittance  function  that  incorporates  the 
averaged  transmittance  is  evaluated  for  each  absorber-pressure- 
temperature  combination  simultaneously  with  a two-parameter  representation 
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of  the  transmittance  function  solely  in  terms  of  the  atmospheric  variables. 
The  latter  representation  is  generalized  to  include  temperature -dependent 
scaling  coefficients  that  cause  all  data  points  to  cluster  along  a single 
curve.  Microwave  transmlttances  are  not  modeled  in  the  same  way;  only 
monochromatic  microwave  transmlttances  are  computed  in  direct  applications 
for  this  study. 

The  third  study  phase  consists  of  the  calculation  of  transmittance 
weighting  function  profiles  and  outgoing  radiances  appropriate  to  actual 
prescribed  atmospheric  structures.  Weighting  functions  are  graphed  for 
interactive  appraisal.  Transmlttances  for  multiple  absorbers  in  a given 
.nterval  are  treated  as  Independent,  so  that  a total  weighting  function 
is  based  on  the  product  of  all  transmlttances.  Radiances  are  compiled  in 
units  commonly  used  for  satellite  data  archiving,  and  are  also  converted 
to  equivalent  Planckian  temperatures  to  enable  more  meaningful  inter- 
channel  comparisons. 

A most  Important  task  is  the  analysis  of  the  effects  on  the  satellite- 
sensed  radiances  of  variations  in  spectral  position,  viewing  angle,  atmo- 
spheric structure,  and  surface  height  and  radiative  temperature.  The  doc- 
umentation of  the  magnitude  of  these  variations  in  radiance  is  useful  for 
analyses  of  background  radiation  for  other  probing  missions  and  for  plan- 
ning or  design  purposes. 

The  final  step  of  the  study  is  an  overview  of  the  numerical  results 
depicting  background  variabilities.  Interpretative  speculations  or  esti- 
mations of  gross  cloud  effects  can  be  made  by  special  analysis  of  tabu- 
lated results. 


IV  ABSORPTION  COEFFICIENTS 


A.  Infrared 

All  of  the  Infrared  absorption  line  parameters  (with  the  exception 
of  line  shape)  required  for  the  detailed  computation  of  absorption  coeffi- 
cients were  taken  from  the  AFGL  compilations  stored  on  magnetic  tape.8 
A program  labeled  RDGOAT  was  written  to  read  the  AFGL  tape,  extract  the 

line  parameters  [strength  (S),  resonant  position  (v  ),  half-width  (v  ) 

o 1 o 

and  energy  of  the  lower  state  (E")]  for  the  particular  absorbers  and  spec- 
tral regions  of  Interest.  Extracted  data  were  printed  and  stored  on 
another  tape.  Ratios  of  Lorentz  to  Doppler  line  widths  were  compiled  at 
two  temperatures  and  included  In  the  listing.  These  ratios  could  be  used 
If  required  in  approximation  methods  for  computing  absorption  coefficients 
for  Voigt  line  shapes  (cf  Kielkopf3),  given  the  results  for  Lorentz  line 
shapes. 

For  this  investigation,  which  was  confined  to  troposphere  influences, 
it  was  considered  adequate  to  use  only  the  Lorentz  line  shape  even  though 
the  weighting  function  for  some  intervals  extended  into  the  stratosphere. 
Thus  the  absorption  coefficient  is  given  by 

k(v)  = Sy/rt[(v  - vq)2  + y2J  . (1) 

The  pressure  (P)  and  temperature  (T)  dependence  of  the  half-width  can  be 
described  by 


Y “ Y0(*/V  (T0/T)"  (2) 

where  the  exponent  n is  usually  taken  as  0.5  or  0.62,  depending  on  the 
absorber  and  type  of  broadening.  The  temperature  dependence  of  the  line 
strength  can  be  approximated  with 

S - (So  Qy(To)  Qr(To)/[Qv(T)  Qr  (T)]}  exp  [1.439  E"(T-Tq)/T  Tq]  (3) 


where  the  vibrational  and  rotational  transition  functions  also  have  been 
approximated  by 


Qv(T)  = a + bT  + cT 


Q (T)  - (T/T  )J  • 


The  exponent  j is  1.5  for  H^,  0^  and  CH^  for  other  absorbers  In  this 
study,  j Is  unity.  The  constant  a Is  assigned  the  value  unity  while  b 
and  c are  taken  as  zero  with  the  exceptions  defined  in  Table  1. 


Table  1 

CONSTANTS  USED  IN  APPROXIMATION  OF 
TEMPERATURE  DEPENDENCE  OF  LINE  STRENGTHS 
(see  Eqs.  4 in  text) 


Molecule 


1.05453 

1.06283 

1.056 


-8. 173E-04  3.20E-06 

-7.422E-04  2.31E-06 

-8. 91 IE -04  3.82E-06 


Given  the  absorber  and  the  wavenumber,  only  the  pressure  and  temperature 
must  now  be  specified  for  a computation  of  absorption  coefficient.  To  be 
compatible  with  files  of  absorption  coefficients  generated  previously, 
computations  were  made  for  three  temperatures  (300,  250,  and  200  K)  and 
four  pressures  (0.9,  0.5,  0.2,  and  0.06  atm  for  ^0;  0.9,  0.3,  0.09,  and 
0.03  atm  for  uniformly  mixed  gases;  0.9,  0.25,  0.08,  and  0.025  atm  for 
0^).  At  each  computational  position,  contributions  from  lines  within 
± 10  cm'1  were  included.  The  spacing  between  wavenumbers  at  which  the 
monochromatic  absorption  coefficients  were  compiled  was  varied  from 
,05 at  the  highest  pressure  to  .004  cm  ^ at  low<v  pressures  (still  with 
Lorentz  line  shape) . 


When  combined  with  previous  absorption  coefficient  computations,  all 
significant  absorption  between  20  cm  ^ and  2560  cm  ^ has  been  covered 
with  one  major  exception:  the  ozone  absorption  bands  between  about  990 
and  1110  cm  . Inasmuch  as  the  principal  ozone  absorption  occurs  at 
altitudes  above  the  clouds,  this  spectral  region,  with  a very  large  num- 
ber of  lines,  was  excluded  for  reasons  of  economy.  Absorption  lines  also 
were  ignored  in  adjacent  regions  dominated  by  the  water  vapor  continuum. 


Absorption  coefficients  were  generated  for  a wide  variety  of  inter- 
vals and  absorbers;  a number  of  spectral  gaps  from  previous  work  were 
filled.  In  addition  to  1^0,  CO^.  and  some  0^  bands,  other  absorbers  in- 
cluded N^O,  CO,  and  CH^.  The  latter  three  gases  were  assumed  to  be  uni- 
formly mixed,  as  was  C02.  Consequently,  it  was  possible  to  condense  the 
absorption  coefficients  in  regions  where  more  than  one  of  the  uniformly 
mixed  gases  was  active.  If  each  absorption  coefficient  is  multiplied 
(weighted)  by  the  ratio  of  its  mixing  ratio  to  that  for  CO^,  then  the 
weighted  sum  of  the  combined  absorption  coefficients  can  be  treated  as 
the  new  absorption  coefficient  for  an  equivalent  concentration  of  CO^. 

If  each  of  the  mixing  ratios  is  denoted  by  r,  then,  for  example, 


where  the  single  combination  absorption  coefficient  (UMG  denotes  Uni- 
formly Mixed  Gas)  is  applied  only  to  equivalent  amounts  of  carbon  dioxide 
along  the  path. 


Figure  2 illustrates  a resultant  UMG  absorption  coefficient  distri- 
bution that  originates  from  a combination  of  C02  and  N20  absorption  co- 
efficients. It  is  apparent  that  the  C02  distribution  is  very  similar  to 
the  final  UMG  distribution,  despite  the  greater  magnitude  of  the  absorp- 
tion coefficients  for  N^.  The  reason  Cor  this  is  that  the  1^0  mixing 
ratio  in  the  atmosphere  is  three  orders  of  magnitude  less  than  that  of 


FIGURE  2 ABSORPTION  COEFFICIENTS  FOR  INDIVIDUAL  GASES  AND  THEIR 
UNIFORMLY  MIXED  ATMOSPHERIC  COMBINATION  (TEMPERATURE 
PRESSURE  - 0.9  atm) 


CO^.  Only  the  stronger  lines  of  N^O  below  about  1890  cm  * have  an  Influ- 
ence on  the  UMG  curve;  easiest  to  detect  are  lines  near  1872,  1880,  and 
1886  cm 

Absorption  coefficients  for  the  water  vapor  continuum  and  the 
nitrogen-induced  continuum  In  this  study  have  been  taken  directly  from 
the  LOWTRAN  3B  model  of  Selby  et  al.4  Extensions  of  coefficients  (for 
both  self  and  foreign  gas  broadening  and  their  temperature -dependent  co- 
efficients) out  to  325  cm  ^ within  the  rotational  band  of  H-0  were  made 

* 1 

possible  by  preliminary  Information.  Adjustments  appearing  In  LOWTRAN 
3B  are  discussed  primarily  by  Roberts  et  al.6 


B.  Microwave 

Whereas  the  monochromatic  absorption  coefficients  computed  for  ab- 
sorbers In  the  Infrared  were  used  In  modeling  over  Intervals  as  large  as 
20  cm  , the  absorption  coefficients  for  the  microwave  were  both  computed 
and  applied  monochromatlcally . Thus,  for  the  microwave  region  no  flies  of 
absorption  coefficients  were  prepared  and  stored;  each  coefficient  Is 
computed  as  needed  in  application. 

The  absorption  coefficient  throughout  the  5 mm  oxygen  band  was  based 
entirely  on  the  treatment  presented  by  Rosenkranz6  with  very  useful  ap- 
proximations for  dealing  with  the  theory  of  overlapping  lines.  The  ab- 
sorption coefficient  k(v)  may  be  expressed  by 

k(v)  = CP(v/T)2F  (6) 

where  C Is  a constant.  The  shape  function  F may  be  described  by 

p ■ plf®»  K(v)  + 4(-v)  + + fN(-v)] + ■7v(''2  + p2“b)}  (7> 

where  |u  is  the  fractional  population  of  the  Initial  state  associated 
with  rotational  quantum  number  N,  and  the  summation  (over  odd  values  of 


* 

Courtesy  of  Dr.  Robert  M.  McClatchey,  AFGL 
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N)  of  the  shape  function  factors  f covers  all  branches  of  significant 
lines  (including  those  at  negative  symmetric  frequencies)  contributing  to 
resonant  absorption.  The  last  term  represents  the  nonresonant  contribu- 
tion, with  wfa  representing  nonresonant  half -width  per  unit  pressure. 

Other  near-diagonal  elements  (w)  of  the  transition  rate  matrix,  frequen- 
cies and  amplitudes  of  oxygen  lines,  and  interference  coefficients  are 
ail  provided  in  tabular  form  by  Rosenkranz;  each  of  the  shape  function 
factors  is  defined  in  terms  of  these  data.  Temperature  dependence  enters 
through  the  fractional  population  and  through  the  resonant  and  non- 
re  sonant  line  half -widths. 

For  water  vapor  at  the  microwave  frequencies,  the  absorption  coeffi- 
cient was  computed  using  two  different  approaches  and  two  line  shapes. 

One  expression  for  the  absorption  coefficient  was  discussed  by  Staelin;7 
the  other  approach  was  described  by  Gaut  and  Reifenstein.8  The  Staelin 
and  Gaut  expressions,  respectively,  for  the  monochromatic  absorption  co- 
efficient (used  here  in  terms  of  nepers  km  *)  can  be  put  in  the  forms 

c(v)  = (pw/'r2*5)v2Av[cie'c/TF  + C2t]  (8a) 

and 

o(v)  = (p^T2*5)*2^^'071?  + C^T1*"]  (8b) 

where  av  - C3  (P/Tn)(l  + C4pwT/P)  and  the  constants  C2>  C3>  and 
differ  with  model  and  resonant  frequency.  F is  the  line  shape  factor, 

F “ \ J + 1 2 

(v-vQ)  + Av  (V+VQ)  + Av 

for  the  Van  Vleck-Weisskopf  shape,  p t8  the  water  vapor  density,  c in- 

w 

eludes  the  energy  of  the  initial  state,  and  av  is  the  appropriate  line 
half-width.  With  the  Gaut  model  and  the  resonant  frequency  of 
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V INFRARED  TRANSMITTANCE  MODELING 


The  computed  monochromatic  absorption  coefficients  for  pertinent 

gases  at  realistic  ranges  of  pressure  and  temperature  form  the  basis  for 

transmittance  modeling.  The  modeling  procedure  is  the  same  as  that  ap- 
1 * 

plied  previously.  Fundamentally,  the  transmittance  model  sought  is  a 
simple  formulation  (involving  a minimum  of  parameterization)  that  would 
be  expressable  explicitly  in  terms  of  atmospheric  variables  and  should 
provide  a sufficiently  universal  fit,  regardless  of  temperature  or  pres- 
sure, that  all  points  follow  a single  curve.  If  a good -enough  fit  is 
achieved,  then  it  is  possible  to  replace  the  nonhomogeneous  path  with  an 
equivalent  homogeneous  path.  Previous  experience  with  the  model  technique 
suggested  that  the  best  results  might  be  obtained  for  spectral  intervals 
of  about  15  cm  or  greater. 

To  begin  the  transmittance  modeling  (see  Figure  3)  spectral  intervals 
must  be  selected  and  a meaningful  set  of  absorber  amounts  must  be  provided 
for  association  with  the  absorption  coefficients.  In  practice,  five  dif- 
ferent absorber  amounts  were  assigned  to  each  of  four  pressures  at  each  of 
three  temperatures.  The  amounts  were  selected  to  conform  to  realistic 
encounters  over  atmospheric  paths.  If  the  absorption  was  too  weak,  the 
amounts  would  be  increased  to  provide  sufficient  data  samples  for  the 
computational  phases. 

The  first  computational  step  is  to  compute  all  of  the  monochromatic 

transmittances  associated  with  absorption  coefficient  k and  absorber 

v 

amount  W: 


T » exp  [ -k  (P,T)  Wl 
v v 


(10) 


Previous  results  could  not  be  used  here  because  they  were  designed  for 
specific  instrument  response  functions. 
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The  many  monochromatic  transmittances , after  convolution  with  a 


spectral  response  function  uj  , are  compressed  by  several  orders  of  magni- 
tude through  the  integration  over  the  given  spectral  interval  av 


T 

Av 


/T  OU 
V V 


(11) 


For  this  study  a square  response  function  of  unit  magnitude  over  the  in- 
terval Av  was  used.  [At  this  point  only  a single  absorber  is  treated  at 
one  time  for  the  interval  Av.  When  final  results  are  applied,  all  sig- 
nificant absorbers  within  a spectral  Interval  are  treated  as  Independent, 
and  the  total  transmittance  is  given  as  the  product  of  the  individual 
transmittances.]  Inasmuch  as  the  unit  response  function  was  used  in  this 
study,  it  was  possible  to  define  the  polychromatic  averaged  transmittances 
within  subintervals  (say  5 cm  S with  the  capability  for  simple  extension 
to  any  larger  interval  that  was  some  multiple  of  the  subinterval.  Accord- 
ingly, an  initial  set  of  intervals  (multiples  of  5 cm  S and  absorbers 
were  identified  (see  Table  2)  for  the  polychromatic  computations.  It  was 
from  subsets  of  these  data  that  polychromatic  transmittances  were  input  to 
the  remaining  modeling  stages. 

As  mentioned  previously,  the  absorption  coefficients  for  all  gases 
considered  to  be  uniformly  mixed  in  the  atmosphere  (e.g.,  CO,,  CO,  CH^, 

NjO)  were  combined  and  treated  as  a single  absorber  (UMG)  with  amounts 
specified  in  terms  of  CO^  amount.  The  UMG  label  included  the  intervals 
with  CO^  only  or  intervals  with  other  uniformly  mixed  gases  in  the  absence 
of  CO^.  Thus,  the  final  infrared  absorbers  with  identified  polychromatic 
transmittances  were  either  H^O,  UMG,  0^,  or  continuum  (discussed  later). 

Given  the  polychromatic  transmittances  for  a variety  of  absorber 
amounts,  pressures,  and  temperatures , the  next  step  was  to  define  a suit- 
able transmittance  function  (i.e.,  some  function  which  depends  only  on  t) . 
Given  the  transmittance  T,  the  transmittance  function  takes  on  a numerical 
value;  on  the  other  hand,  given  the  value  of  the  transmittance  function, 
the  transmittance  can  be  recovered.  An  example  of  a transmittance  function 

might  be  log  t,  or  log  (-log  t) . Such  simple  transmittance  functions  sel- 
e e e 

20 


i 


Table  2 


SPECTRAL  INTERVALS  AND  ABSORBERS  CONSIDERED  IN 
POLYCHROMATIC  TRANSMITTANCE  COMPUTATIONS 


Spectral  Range 

-1 
cm 


Absorbers 


530 

- 

560 

h2o,  (co3),  (n2o) 

680 

- 

7 20 

h2o,  co2,  o3 

720 

- 

760 

h2o,  co2,  o3 

880 

- 

900 

Continuum 

1190 

- 

1210 

Continuum 

* 

1240 

- 

1270 

h2o,  n2o,  ch4 

1360 

JL. 

- 

1390 

h2o,  (ch4) 

1550 

- 

1570 

H2° 

1680 

- 

1700 

h2° 

1920 

- 

1950 

H20,  <C02) 

*2000 

- 

2030 

h2o,  (C02),  (CO) 

2060 

- 

2090 

H20,  C02,  03 

2090 

- 

2120 

h2o,  co2,  03 

2170 

- 

2200 

h2°.  C02.  N20,  CO 

2230 

- 

2260 

h2o,  co2,  n2o,  CO 

*2280 

- 

2310 

(h20),  co2,  (n20) 

2500 

- 

2520 

Continuum 

"ft 

Subsequently  dropped  from  further  analysis  to  reduce  scope. 

^Absorbers  in  parentheses  not  actually  used  in  modeling.  Continuum 
considered  over  entire  spectral  range. 


dom  would  apply  to  all  ranges  of  absorber-pressure-temperature  combinations 
without  introduction  of  a series  of  parameters.  The  form  adopted  here  for 
transmittance  function  is  the  same  as  that  used  and  discussed  previously: 


2 


(12) 


Y 

c 


This  two-term  form  with  the  parameter  C,  where  $ ^ is  the  inverse  error 
function  with  argument  (1  - t) , is  characteristic  of  a combination  of 
regular  and  random  transmission  models,  blended  in  accordance  with  the 
magnitude  of  C.  With  C = 0 the  random  model  prevails  whereas  with  C = 1 
the  regular  model  prevails.  For  almost  all  situations  encountered,  the 
value  of  C ranges  between  1.0  and  -1.0,  but  occasionally  extends  to  -2.0 
with  a good  representation.  When  the  best  fitting  C was  less  than  -2.0, 
the  sensitivity  to  C was  very  slight  or  the  transmittances  did  not  fit 
well  to  a single  curve.  Although  not  encountered  in  practice,  a limita- 
tion of  Eq.  (12)  exists  on  physical  grounds  since  i|rc  must  be  positive. 

For  very  small  T this  function  could  become  negative  for  C greater  than 
unity;  if  five  places  are  retained  in  t and  if  t = .00001,  then  C cannot 
exceed  1.39  for  to  remain  positive.  Even  that  occurrence  could  be  pre- 
vented by  the  addition  of  a suitable  correction  term  including  C,  1 - C, 

2 -1 

and  the  difference  between  in  t and  — { (1-t)  as  factors  (not  required 

in  the  present  application).  ^ 


The  numerical  technique  for  recovering  the  transmittance  T from  the 

transmittance  function  (Eq.  12)  is  given  in  Appendix  B.  What  has  not 

been  discussed  is  the  manner  in  which  the  parameter  C is  determined.  A 

key  step  in  the  modeling  procedure  is  to  solve  for  C simultaneously  with 

the  parameters  of  a representation  of  the  transmittance  function  solely 

★ 

in  terms  of  atmospheric  variables.  The  form  adopted  for  the  atmospheric 
representation,  ^ 

*c  * V / V(W/P)  ^BT2  (13) 


In  a previous  study  the  right  hand  side  of  Eq.  (13)  was  slightly  more 
complicated  with  an  additional  appearance  of  Bt  in  the  denominator 

♦ c “ \W  / ( \Aw/P)  + + Bt  ) . 

For  00  this  form  corresponds  with  the  Malkmas9  model  whereas  Eq.  (13) 
corresponds  with  the  Goody10  model. 

r 
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contains  two  temperature -dependent  parameters,  A^  and  B^,.  In  practice, 
with  C expressed  according  to  Eq.  12,  transmittances  are  introduced  for 
the  reference  temperature  (250  K).  The  equations  are  solved  by  a non- 
linear least-squares  technique  to  determine  the  best  fitting  A^,  , and 

C in  terms  of  a dependent  and  two  independent  variables.  Table  3 contains 
a listing  of  these  constants  for  the  noncontinuum  regions  of  Table  2. 

The  value  of  C so  determined  is  taken  as  a constant  independent  of  tem- 
perature (an  approximation).  Now  the  left  side  of  Eq.  (13)  is  fully  de- 
fined, so  that  the  nonlinear  least-squares  curve-fitting  technique  can  be 
reapplied  to  obtain  the  best  fitting  parameters  A^,  and  B^  at  each  of  the 
three  temperatures  considered  in  the  modeling.  For  some  regions  the  best 
fit  is  obtained  with  the  single  parameter  A^  (i.e.,  let  B^  = 0): 

*c  = AT/WP  . (14) 

The  solution  of  A^,  and  B from  Eq.  (13)  or  A^,  from  Eq.  (14)  has  re- 
sulted in  excellent  representations  of  data  points  at  all  pressures  and 
absorber  amounts  by  a single  curve  for  a given  temperature.  However, 
with  C held  fixed  on  the  left  side,  the  parameters  A^,  and  B^  for  the 
three  temperatures  can  be  used  to  define  temperature -dependent  scaling 
coefficients  for  the  independent  variables,  in  terms  of  a second -degree 
polynomial  in  temperature.  Suppose  the  scaled  independent  variables  in 
Eqs.  (13)  and  (14)  are  defined  as 


* 

(W  ) = 4 

fsL(T)  dW 

(15a) 

(WP)*  = > 

f s2(T)  p dW 

(15b) 

(W/P)*  » 

(W*)2  / (WP)* 

= [/ S3(T)/P  dw]  + 

(15c) 

Not  required  as  separate 


specification  since  S^(T) 


[s^DT 


/ 


S2(T) 
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Table  3 

INITIAL  PARAMETERS  AT  DEFINED  REFERENCE  TEMPERATURE  FOR 
TRANSMITTANCE  MODELING  IN  SELECTED  20  cm-1  INFRARED  INTERVALS 


Interval 


(cm'1) 

*0 

*0 

C 

Water  Vapor 

440-460 

"064 

• 

• a* 

0.199 

446-466 

2.020 

mm 

0.614 

460-470 

1.992 

mm 

0.656 

660-660 

0.762 

mm 

0.747 

656-655 

0.912 

mm 

0.746 

640-660 

0.961 

0.267 

1560-1680 

5.088 

0.817 

1565-1586 

2.405 

— 

0.826 

1570-1590 

2.866 

— 

0.697 

1920-1940 

0.906 

— 

0.607 

1926-1946 

0.894 

mm 

0.626 

1950-1960 

1.248 

— 

0.275 

680-700 

0.554 

-0.162 

0.326 

686-706 

0.564 

0.184 

0.726 

690-710 

0.569 

0.199 

0.662 

696-716 

0.558 

0.256 

0.378 

700-720 

0.222 

-0.266 

0.176 

TaO-756 

0,206 

-0.206 

-1.469 

710-750 

0.118 

0.506 

•1.161 

716-756 

0.100 

0.460 

-0.587 

720-740 

0.095 

-0.419 

-0.364 

726-746 

0.176 

-0.574 

-0.550 

750-760 

0.162 

-0.364 

0.442 

756-766 

0.177 

-0.384 

0.286 

740-760 

0.181 

-0.587 

0.515 

2060-2080 

0.584 

0.190 

-0.043 

2065-2086 

0.219 

0.181 

0.644 

2070-2090 

0.198 

0.596 

1.029 

2090-2110 

0.162 

0.337 

0.517 

2096-2116 

0.166 

0.652 

0.465 

2100-2120 

0.167 

0.440 

0.180 

2170-2190 

0.077 

1.142 

0.461 

2176-2196 

0.066 

1.39? 

0.616 

2180-2200 

0.064 

1.377 

0.391 

2250-2260 

0.066 

5.420 

-0.005 

2256-2266 

0.066 

5.520 

0.543 

2240-2260 

0.095 

2.614 

-0.922 

Interval 


(c.-1) 

*0 

B0 

C 

Uniformly  Mixed  Oaaea 

680-700 

0.068 

0.623 

0.662 

686-706 

0.622 

0.644 

0.528 

690-710 

0.380 

0.860 

0.502 

696-716 

0.268 

1.256 

0.631 

700-720 

0.228 

2.076 

0.494 

720-740 

0.164 

2.647 

-2.060 

726-746 

0.096 

3.481 

0.779 

730-760 

0.091 

3.746 

0.534 

736-766 

0.070 

4.446 

0.642 

740-760 

0.059 

6.113 

-0.089 

2060-2060 

0.031 

10.444 

0.0 

2066-2086 

0.028 

10.983 

0.0 

2070-2090 

0.023 

10.643 

0.0 

2170-2190 

0.044 

17.624 

0.616 

2175-2196 

0.065 

11.766 

0.454 

2180-2200 

0.070 

9.523 

0.389 

2230-2260 

0.240 

4.230 

0.679 

2236-2266 

0.235 

3.200 

0.792 

2240-2260 

0.230 

1.906 

0.886 

Oaone 


680-700 

0.964 

5.764 

-0.090 

686-705 

0.868 

5.985 

-0.299 

690-710 

0.767 

5.664 

-0.238 

695-715 

0.834 

6.846 

-0.190 

700-720 

0.966 

6.499 

-0.260 

720-740 

1.081 

6.864 

-0.539 

726-746 

0.971 

7.578 

-0.364 

750-760 

0.786 

7.394 

-0.410 

736-765 

0.742 

7.485 

-0.467 

740-760 

0.664 

7.666 

-0.307 

2060-2080 

0.502 

5.762 

-0.133 

2066-2066 

0.804 

3.638 

-0.329 

2070-2090 

1.066 

3.416 

-0.167 

2090-2110 

1.860 

3.006 

-0.278 

2096-2116 

1.850 

5.221 

0.490 

2100-2120 

2.067 

3.298 

-0.637 
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where 


S^T)  = (AtBq)  / (AqBt)  2:  1 + a^T  - Tq)  + a^T-T^2  (16a) 

S2(T)  - A2  / A2  2=  1 + a3(T  - To)  + a4(T  - Tq)2  (16h) 

| S 3 (T ) - B2  / B2  i 1 + a5(T  - Tq)  + a6(T  - Tq)2]+  . (16c) 

The  coefficients  In  Eq.  (16)  are  determined  from  Information  at  the  two 
nonreference  temperatures.  When  the  scaling  coefficients  are  applied  to 
the  variables,  all  data  points  should  fall  on  a single  curve,  rather  than 
a separate  curve  for  each  temperature.  When  this  is  accomplished  a "uni- 
versal" representation  has  been  achieved,  in  which  the  single  curve, 

iH  ! AW*  / \ (W/P)*  + B2  (17) 

describes  the  transmittance  regardless  of  absorber  amount,  pressure,  or 
temperature.  If  such  a universal  curve  has  been  achieved,  then  it  has 
been  demonstrated  that  a corresponding  equivalent  homogeneous  path  can 
replace  the  nonhomogeneous  path.  The  explicit  relationship  to  arbitrary 
distributions  of  atmospheric  variables  has  great  potential  for  generating 
matrices  of  coefficients  that  could  be  used  in  fast  techniques1 1*  1 8 * 13 
for  extending  transmittancen  through  the  atmosphere. 

Although  success  has  been  achieved  with  the  modeling  technique  de- 
scribed, that  is  not  to  say  that  the  technique  will  always  work  anywhere 
or  that  adequate  empirical  tests  have  been  made.  In  fact,  for  some  spec- 
tral regions  (not  considered  here)  problems  were  experienced.  A number  of 
alternative  forma  to  that  expressed  hy  Eq.  (13)  or  (14)  were  tried,  but 
generally  met  with  less  success.  The  lack  of  improvement  persisted  even 
with  the  addition  of  a single  term  and  another  parameter. 


See  footnote  on  previous  page. 


25 


Examples  of  the  fits  achieved  by  Eqs.  (13)  and  (14)  reveal  how  well 
the  theoretical  data  can  be  modeled.  In  the  illustrations,  different 
symbols  are  maintained  for  the  three  different  temperatures,  but  the  var- 
ious pressures  and  absorber  amounts  within  each  temperature  category  are 
not  distinguished.  Figure  4 illustrates  the  results  obtained  from  Equa- 
tion (14)  for  water  vapor  in  the  interval  530-550  cm.  At  the  top,  data 
points  are  plotted  as  computed,  but  the  curve  defined  for  standard  tem- 
perature is  shifted  in  accordance  with  the  scaling  factors  for  300  K and 
200  K.  Each  shifted  curve  and  the  reference  curve  (250  K)  follow  the  dis- 
tribution of  points  very  closely.  Similarly,  when  the  universal  curve  at 
the  reference  temperature  is  plotted  and  data  points  are  shifted  in  ac- 
cordance with  Equations  (15)  and  (16),  then  a single  curve  with  well 
clustered  points  is  obtained,  as  shown  in  Figure  4b. 

Similar  results  are  obtained  for  water  vapor  in  another  spectral  in- 
terval, 725-745  cm  1 , where  the  absorption  is  much  weaker.  In  this  re- 
gion the  two-parameter  form  of  Equation  (13)  is  superior.  Figure  5a  shows 
the  excellent  results  from  shifting  the  curves  to  meet  the  points  and 
Figure  5b  from  scaling  the  points  to  coincide  with  the  universal  curve. 

It  would  be  difficult  to  improve  the  fit.  Figure  6a  shows  the  fit  of  all 
data  points  (with  temperature-dependent  shifts  applied)  to  the  reference 
curve  for  ozone  in  the  same  spectral  region.  Here  the  data  points  occur 
in  distinct  clumps  (determined  by  choice  of  conditions)  as  a result  of  the 
lack  of  any  significant  shifts  with  temperature,  which  normally  would 
shift  warmer  points  in  one  direction  and  colder  points  in  the  opposite  di- 
rection. Figure  6b  shows  a similar  plot  for  CO^  for  which  more  dispersion 
results  from  the  temperature -dependent  scaling.  It  may  be  noted  that  some 
points  on  one  tail  of  the  curve  do  not  fit  with  the  curve.  This  represents 
about  the  most  serious  departure  observed  in  any  of  the  fitted  curves. 

The  slight  displacement  probably  is  due  to  the  Q branch  near  741  cm'1 
A plot  of  the  C02  absorption' coefficients  for  the  spectral  region  is 
shown  in  Figure  7. 
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FIGURE  6 UNIVERSAL  REPRESENTATIONS  FOR  OZONE  AND  CARBON  DIOXIDE 

TRANSMITTANCES.  726-746  cm*1  (FOUR  PRESSURES,  THREE  TEMPERATURES) 
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An  initial  fit  for  the  spectral  region  720-740  cm  ^ led  to  a satis- 
factory fit  but  the  resultant  parameters  were  quite  different  from  those 
of  neighboring  intervals.  The  difficulty  was  isolated  to  the  720-725 
portion  of  the  interval.  That  portion  of  the  interval  is  significantly 
influenced  by  the  720  cm  Q branch,  while  the  absorption  coefficients 
over  the  rest  of  the  interval  show  a distinctly  different  distribution. 

It  is  believed  that  since  the  absorption  in  the  two  regions  arises  from 
different  portions  of  the  atmosphere,  it  is  difficult  to  portray  atmo- 
spheric influence  in  the  response  of  a single  channel.  One  way  to  avoid 
this  type  of  problem  would  be  to  use  a separate  narrow  channel  to  cover 
the  Q branch  and  another  channel  (or  more  narrow  channels14)  to  cover  the 
remainder.  Thus,  it  would  be  best  to  select  the  intervals  after  examin- 
ing plots  of  all  pertinent  absorption  coefficients  or  of  computed  specific 
spectra.16  In  this  study  of  20  cm  intervals,  attention  was  merely 
shifted  from  720-740  cm  ^ to  725-745  cm  Although  a very  good  fit  was 
achieved  with  the  shift,  it  can  be  seen  that  the  weaker  Q -branch  was  picked 
up  at  the  expense  of  the  strong  one. 


Results  of  the  modeling  for  spectral  intervals  retained  in  this  study 
are  summarized  in  Table  4.  Included  are  the  constant  C of  the  transmit- 
tance function,  the  parameters  AQ  and  BQ  of  the  atmospheric  representation, 
and  the  temperature -dependent  scaling  coefficients  of  Eq.  (17). 


All  of  the  transmittance  modeling  discussed  to  this  point  has  dealt 
only  with  absorption  lines--not  the  continuum.  Although  the  method  for 
handling  the  continuum  is  described  by  Selby  et  al.,4  the  general  expres- 
sions for  the  H20  continuum  can  be  put  in  the  form 


T ■ 


(18) 


where  kg  is  the  absorption  coefficient  at  temperature  296K.  for  self- 
broadening  and 


U-O.lp  !pw  exp[C2(296/T)  -1)]  + C3(P  - Pw)  exp[C4(296/T)  - 1)]J 

W l W W (19) 
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for  n (water  vapor  density)  in  g m , p (vapor  pressure)  and  P (total 
r w w 

pressure)  in  atm,  and  T (temperature)  in  Kelvins.  The  height  z is  taken 
in  km.  In  application,  the  constant  either  takes  on  the  value  0 or 
C2>  but  it  can  take  on  other  independent  values.  For  the  nitrogen  con- 
tinuum, the  transmittance  has  been  expressed  as 


In  previous  treatments  of  the  continuum  we  have  used  exponential  functions 
similar  to  those  in  Eqs.  (18)  and  (20)  but  with  the  argument  raised  to 
some  power  less  than  unity. 


Table  4 

P I 

SELECTED  PARAMETERS  FOR  UNIVERSAL  TRANSMITTANCE  MODEL 
AND  POR  CPUTINUUM  REPRESENTATION 


Interval 

(cm-1) 

C 

A 

B2 

*1 

*2 

*3 

445-466 

0.5 

2.03134 

Water  Vapor 

0.0 

2.174 

1.581 

530-550 

0.7 

0.79276 

0.0 

— 

- — 

1.362 

0.672 

700-720 

0.2 

0.22150 

0.06542 

2.518 

2.081 

2.471 

1.935 

720-740 

-0.4 

0.09310 

0.17687 

1.802 

0.973 

2.442 

2.151 

725-745 

-0.4 

0.17643 

0.14052 

2.756 

2.376 

2.814 

2.531 

735-766 

0.3 

0.17725 

0.14711 

3.263 

3.207 

3.053 

2.883 

1365-1385 

0.8 

2.47774 

0.0 

• • • 

... 

1.051 

0.186 

1925-1945 

0.6 

0.90899 

0.0 

- — 

mmm 

1.052 

0.224 

2065-2065 

0.8 

0.21990 

0.03344 

1.215 

-0.069 

1.817 

1.003 

2180-2200 

0.4 

0.06421 

1.89354 

2.524 

2.019 

3.314 

3.328 

2240-2260 

-0.9 

0.09309 

6.29549 

2.770 

2.377 

3.695 

4.067 

700-720 

0.5 

Uniformly  Mixed  Gases 

0.22782  4.29734  1.841 

0.995 

2.433 

1.876 

720-740 

-2.0 

0.16327 

6.43306 

1.521 

0.476 

1.946 

1.097 

725-745 

0.8 

0.09507 

12.05173 

1.890 

1.044 

2.210 

1.689 

735-756 

0.5 

0.07078 

19.94220 

2.351 

1.698 

2.883 

2.681 

2065-2065 

0.0 

0.02747 

120.62780 

-0.058 

0.043 

0.328 

0.229 

2180-2200 

0.4 

0.06970 

90.38351 

0.617 

0.025 

0.937 

-0.162 

2240-2260 

0.9 

0.22877 

3.60314 

0.933 

0.141 

1.244 

0.634 

700-720 

-0.3 

0.96834 

Ozone 

30.99536 

-0.089 

-0.027 

-0.164 

-0.019 

720-740 

-0.6 

1.06900 

46.08662 

0.122 

-0.095 

0.167 

-0.114 

726-745 

-0.4 

0.98043 

58.34086 

0.134 

-0.062 

0.260 

-0.129 

736-766 

-0.5 

0.74925 

66.89678 

0.114 

-0.061 

0.227 

-0.116 

2065-2085 

-0.3 

0.79759 

14.48916 

0.416 

-0.167 

0.803 

-0.227 

Continuum 


446-465 

166.98 

1.33 

0.014 

-3.72 

... 

630-660 

71.81 

3.60 

0.007 

-2.16 

... 

700-720 

23.881 

5.60 

0.002 

0.0 

aee 

T80-740 

21.643 

6.66 

0.002 

0.0 

726-746 

20,986 

5.68 

0.002 

0.0 

... 

736-765 

19.673 

6.71 

0.002 

0.0 

... 

680-900 

9.186 

6.08 

0.002 

0.0 

ammmm 

1190-1210 

4.609 

6.08 

0.002 

0.0 

... 

1366-1385 

4.176 

6.08 

0.002 

0.0 

... 

1926-1946 

0,67 

4.66 

0.12 

0.0 

... 

2066-2066 

0.64 

4.66 

0.12 

4.66 

.00024 

2180-2200 

0.40 

4.66 

0.12 

4.66 

.02890 

2240-2260 

0.33 

4.66 

0.12 

4.66 

.06595 

2600-2620 

0.118 

4.66 

0.12 

4.66 

.02469 

33 


fA* 

{ J 


VI  WEIGHTING  FUNCTIONS 


An  objective  of  the  study  was  to  examine  several  spectral  intervals, 
distinctly  separated  in  wavenumber,  with  similar  transmittance  weighting 
to  the  tropospheric  regions  where  clouds  are  common  and  to  the  lower 
boundary.  Differences  (or  lack  of  differences)  in  the  multispectral  radi- 
ances might  then  provide  background  information  that  would  be  useful  in 
subsequent  remote  sensing  of  clouds  themselves. 

The  specification  of  widely  separated  spectral  intervals  having 
identical  transmittance  weighting  factors  for  the  viewed  earth/atmosphere 
is  difficult  to  achieve,  at  least  without  the  introduction  of  deliberate 
filter  designs  for  the  receiver.  Without  a more  sophisticated  system  ap- 
proach than  that  possible  here  it  is  difficult  to  find  a precise  match  in 
response  for  two  different  intervals.  Unless  the  absorbers  played  the 
same  relative  role  within  the  two  intervals,  it  would  be  likely  that  if  a 
match  was  obtained  for  one  condition  (atmospheric  structure  or  angle  of 
view)  the  match  would  not  hold  when  that  condition  changed.  Furthermore, 
even  if  general  similarity  in  attenuation  occurred  as  a result  of  distri- 
butions of  absorption  lines,  there  would  still  be  a requirement  that  the 
background  continuum  in  each  interval  be  the  same  (at  least  for  low-level 
responses),  which  is  relatively  uncommon.  At  this  stage  of  study,  for 
intervals  of  about  20  cm  ^ , only  approximate  matches  in  response  were  re- 
quired to  obtain  pertinent  background  information. 

An  obvious  first  step  in  the  selection  of  possible  matching  intervals 
in  different  parts  of  the  spectrum  ir  to  examine  the  absorption  coefficient 
distributions  (cf.  Figures  2 and  7),  but  this  is  difficult  since  there  can 
be  three  different  types  of  absorbers  present,  each  with  varying  atmo- 
spheric distributions.  Selections  also  can  be  based  on  the  parameters 
derived  in  the  transmittance  modeling  (cf.  Table  3),  but  again  each  set 
of  parameters  refers  separately  to  each  type  of  absorber  and  does  not  in- 
clude the  continuum.  Both  approaches  were  used  in  this  study,  with  the 
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limitation  that  matching  be  done  for  infrared  intervals  as  wide  as  20  cm 
or  more.  Approximate  matches  could  be  made  more  nearly  identical  by 
introducing  an  artificial  spectral  filter  that  would  force  the  match  in 
response  (whether  or  not  such  a filter  actually  could  be  constructed). 
Useful  information  could  be  derived  in  this  way,  neglecting  instrument 
design,  which  is  a later  step. 

A final  check  on  the  match  between  two  intervals  is  accomplished  by 
comparing  the  transmittance  weighting  functions  themselves.  This  check- 
ing is  done  whether  or  not  any  improvement  can  be  made.  To  be  useful, 
the  check  must  be  made  with  the  total  transmittance  weighting  function 
for  different  conditions.  The  total  transmittance  is  defined  by  the  prod- 
uct of  the  separate  transmittances  for  each  of  N absorbers  and  the  trans- 
mittance for  the  continuum;  that  is, 

N+l 

T = 7 T T (i)  . (21) 

i-1  Av 

To  appreciate  the  choice  of  weighting  function  it  is  helpful  to  de- 
scribe the  radiance  R(v^)  for  the  ith  spectral  interval  as  measured  from 
above  the  atmosphere  at  height  H 


R( 


V ■ B<vV  T<vs)  -/.»S  B<vV  ^ 


dz 


or,  by  dropping  the  spectral  identification  and  rewriting, 


R - B(S)t(S)  +J  B(z)  dz 


(22) 


where  the  total  transmittance  T is  always  measured  over  the  path  from  H 
to  the  level  indicated  by  the  argument.  The  Planck  function  B(z)  is  de- 
fined in  terms  of  wave  number  v by 


B(z)  - C1vJ/(exp[C2v1/T(z)]  - l) 

where  C,  and  C are  constants. 
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(23) 
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For  this  study  dT(z)/dz  has  been  selected  as  the  transmittance 
weighting  function,  with  units  km  The  first  term  on  the  right  side  of 
Eq.  (22)  depicts  the  contribution  from  the  surface  emission  (assumed  here 
to  be  blackbody)  that  is  transmitted  through  the  entire  atmosphere.  The 
integral  represents  the  contribution  from  the  molecular  atmospheric  emis- 
sion, after  being  transmitted  from  height  z to  level  H.  Actual  radiance 
contributions  depend  on  the  Planck  function  and  the  temperature  distribu- 
tion; therefore,  profiles  of  the  transmittance  weighting  function  describe, 
for  the  molecular  atmosphere,  the  relative  weighting  of  Planckian  contri- 
butions to  the  observed  radiance. 

The  transmittance  weighting  function  is  the  derivative  of  the  trans- 
mittance with  respect  to  the  height  variable.  Thus,  inspection  of  the 
slope  of  the  transmittance  profile  can  be  interpreted  in  terms  of  the 
weighting  function.  Figure  8 presents  the  transmittance  and  the  weight- 
ing function  profiles  at  relatively  weak  and  strong  absorption  frequen- 
cies of  oxygen  in  the  microwave  region.  At  52.85  GHz  the  trnsmittance 
profile  shows  an  increasing  slope  below  about  20  km  and  a slowing  of  the 
rate  of  decrease  of  transmittance  near  the  ground.  This  translates  to 
the  weighting  function  profile  shown  in  Figure  8,  which  has  a maximum 
near  the  ground.  At  54.9  GHz  the  transmittance  essentially  goes  to  zero 
at  the  ground  surface,  and  the  profile  shows  the  steepest  slope  of  t 
with  respect  to  z at  some  distance  above  10  km.  This  profile  leads  to 
a weighting  function  profile  with  a maximum  just  above  10  km,  with  tails 
approaching  zero  both  at  the  top  and  bottom.  Similarly,  if  the  attenu- 
ation is  increased,  as  a result  of  either  slant  path  viewing  or  an  in- 
crease in  the  absorber,  the  weighting  function  profile  is  raised  somewhat 
in  altitude,  with  some  increase  in  contribution  at  higher  altitudes  and, 
for  opaque  regions,  a decrease  in  contribution  at  the  bottom.  The  char- 
acteristics will  be  evident  regardless  of  the  wavenumber  or  frequency. 

(In  keeping  with  conventional  usage,  we  have  retained  the  frequency  no- 
tation for  the  microwave  regions.) 
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Note  the  smooth  broad  nature  of  the  0?  weighting  function  profiles 
in  Figure  8,  despite  the  fact  Chat  the  computations  are  for  a single  fre- 
quency. The  monochromatic  treatment  of  water  vapor  in  the  microwave  re- 
gion shows  more  irregularity,  in  large  part  because  of  the  irregularity 
of  the  water  vapor  distribution. 

Looking  again  at  the  interval  725-745  cm  \ we  see  that  Figure  9a 
shows  the  transmittance  weighting  functions  for  each  of  the  four  component 
transmittances  as  well  as  for  the  total  for  the  Midlatitude  Summer  atmo- 
spheric profile  of  AFGL.16  Figure  9b  shows  the  same  thing  for  the  Midlati- 
tude Winter  profile.  Perhaps  the  most  impressive  feature  is  the  similarity 
of  the  profile  for  each  total  with  that  for  CO^,  even  though  some  of  the 
other  (individual)  profiles  appear  to  be  significant.  It  is  easy  to  be 
somewhat  misled  by  the  individual  component  profiles,  since  the  total 
weighting  function  is  not  a combination  of  individual  weighting  functions. 
Instead,  the  total  transmittance  must  be  formed,  Eq,  (21),  before  the  de- 
rivative with  respect  to  height  is  determined. 

It  is  obvious  from  Figure  9 that  the  ozone  plays  only  a very  minor 
role,  restricted  to  the  stratosphere.  Both  the  water  vapor  line  absorp- 
tion and  the  continuum  have  much  less  weight  in  winter  than  in  summer  be- 
cause of  the  drier  profile.  Figure  10  shows  the  total  transmittance 
weighting  profile  for  0°  zenith  angles  in  both  summer  and  winter  and  for 
40°  zenith  angle  in  summer.  This  plot  illustrates  variability  in  the 
weighting  function,  resulting  from  changes  in  profile  and  in  viewing  angle. 
Of  course  the  variability  depends  on  the  amount  of  absorber  present  as  well 
as  the  degree  of  attenuation.  Figure  11a  shows  the  variabilities  for  53.5 
GHz  oxygen  attenuation  for  the  same  conditions  shown  in  Figure  10.  Figure 
lib  illustrates  corresponding  variability  for  the  170 -GHz  water  vapor  trans- 
mittance weighting  profile.  Variations  in  the  curve  of  Figure  lib  are 
more  striking  because  of  the  real  seasonal  change  in  water  vapor  and  be- 
cause of  the  strong  low-level  attenuation. 

Figure  12  illustrates,  for  both  summer  and  winter,  the  good  match 
found  for  the  total  weighting  functions  of  the  intervals  445-465  cm  and 
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WEIGHTING  FUNCTION  — km'1 


FIGURE  10  TOTAL  TRANSMITTANCE  WEIGHTING  FUNCTIONS,  728-745  cm 


FIGURE  11  TRANSMITTANCE  WEIGHTING  FUNCTIONS  FOR  53.5  GHz  (OXYGEN)  AND  170  GHz  (WATER  VAPOR) 


HEIGHT  — km  HEIGHT 


FIGURE  12  SUMMER  AND  WINTER  PROFILES  OF  WEIGHTING  FUNCTIONS  FOR 
SIMILAR  SPECTRAL  REGIONS,  466-486  cur1  AND  1366-1385  cm-1 
VERTICAL  PATHS 


1365-1385  cm  , despite  the  much  weaker  continuum  over  the  latter  inter- 
val. These  intervals  respond  to  midtropospheric  water  vapor  with  essen- 
tially no  surface  contribution  in  summer. 

For  less  strongly  absorbing  water  vapor,  Figure  13  shows  an  approxi- 
mate match  between  two  other  20  cm  intervals:  530-550  cm  and  1925-1945 
cm  *■.  The  response  from  530-550  cm  * is  stronger  than  that  for  1925-1945 
cm  , but  the  profile  is  very  similar.  In  this  case,  it  should  be  rela- 
tively easy  to  obtain  an  excellent  match  with  differential  filtering  be- 
tween radiometric  channels. 

The  effect  of  changing  the  bandwidth  from  20  cm  1 to  10  cm  * is 
shown  in  Figures  14a  and  14b  for  both  spectral  regions  of  Figure  13.  It  is 
seen  that  in  both  regions  as  the  10  cm  1 interval  is  shifted  toward  higher 
wavenumber  the  absorption  increases.  The  10  cm  1 profiles  for  1925-1945 
show  blunted  and  relatively  subdued  maxima.  In  both  cases  the  last  in- 
terval differs  in  shape  from  the  first  two  10  cm  1 Intervals.  (Therefore, 
significant  changes  in  the  parameterization  can  result  from  changes  in 
bandwidth. ) 

Better  spectral  matches  have  been  obtained  for  other  Intervals. 

Figure  15  shows  the  good  match  between  530-550  cm  and  the  178  GHz  micro- 
wave  frequency  (water  vapor)  for  the  Midlatitude  Summer  profile.  Figure 
16  shows  the  apparent  match  between  720-740  cm  1 and  2180-2200  cm  * at 
separate  spectral  regions  comnonly  applied  to  remote  sensing.  As  can  be 
seen  later,  the  larger  return  from  the  upper  atmosphere  in  the  720-740 
cm  1 interval  prevents  an  actual  match.  In  fact,  despite  differences  in 
the  water  vapor  contribution,  the  2180-2200  cm  1 interval  matches  better 
with  the  735-755  cm  1 interval,  but  that  latter  Interval  still  yields 
lower  equivalent  temperatures.  Figure  17  shows,  for  both  summer  and 
winter,  the  good  match  between  weighting  function  profiles  for  700-720 
cm”  , 2240-2260  cm"1,  and  54  GHz  (oxygen).  With  only  slight  adjustments 
in  wavenumber  or  frequency  an  even  better  match  could  be  obtained  (the 
oxygen  absorption  is  too  weak  at  54  GHz  to  provide  more  response  to  higher 
altitudes) . 
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FIGURE  13  TRANSMITTANCE  WEIGHTING  FUNCTION  COMPARISON,  630-560  crrr1  AND 
1925-1946  cm-1,  WINTER  PROFILES 


OVERLAPPING  10  cm’1  INTERVALS 


Figure  18  illustrates  transmittance  weighting  profiles  for  three  dif- 
ferent window  regions  for  which  only  the  continuum  was  considered.  For 
the  summer  profile  the  880-900  cm  absorption  exceeds  that  for  1190-1210 
cm  and  2500-2520  cm  . The  latter  region  shows  the  weakest  response, 
but  with  an  influence  over  considerably  more  depth.  This  results  from 
the  nitrogen  continuum  peculiar  to  that  interval.  In  winter  the  2500- 
2520  cm  * region  shows  little  change  while  the  weighting  functions  for 
the  other  two  intervals  are  greatly  reduced.  T.n  fact,  2500-2520  cm 
shows  the  strongest  absorption  for  the  winter  profile. 

Figure  19  illustrates  the  transmittance  weighting  at  three  microwave 
frequencies.  The  solid  curves  show  the  computed  joint  transmittance  for 
H^O  and  0^,  whereas  the  dashed  curves  show  the  curves  for  H^O  only. 
Although  the  attenuation  is  weak  in  all  cases,  the  frequency  22.235  GHz 
corresponds  to  a weak  water  vapor  line.  For  H^O  only,  the  frequency 
37  GHz  shows  a slightly  smaller  weighting  function  than  that  for  19.35 
GHz.  When  oxygen  attenuation  (from  the  60  GHz  band)  is  introduced,  the 
total  weighting  function  is  smallest  at  19.35  GHz,  but  approximately 
twice  as  large  as  for  H^O  alone.  Additional  weighting  functions  for  H^O 
and  0^  microwave  frequencies  are  shown  in  Appendix  C.  The  versatility 
offered  by  the  microwave  region  for  selection  of  weighting  function  pro- 
files is  evident. 


HEIGHT 


VII  RADIANCES 


A.  Upwelling  Radiances 

The  expression  of  Equation  22  for  the  radiance  received  from  above 
the  atmosphere  is  strongly  dependent  on  the  Planck  function.  Because  the 
Planck  function  shows  a marked  variation  in  wavenumber,  especially  for 
temperature  changes,  the  radiances  will  vary  considerably  with  wavenumber. 
At  300  K the  blackbody  maximum  occurs  at  about  600  cm  ^ , whereas  at  200  K 
the  blackbody  maximum  occurs  at  400  cm  In  addition,  at  the  higher 
wavenumbers  the  Planck  function  is  much  more  sensitive  to  temperature 
differences  than  at  low  wavenumbers. 

At  sufficiently  low  frequencies,  the  Planck  function  becomes  a lin- 
ear function  of  temperature.  Consequently,  the  brightness  temperature  of 
the  microwave  radiometer  can  be  substituted  for  the  radiance,  and  the 
Planck  functions  appearing  singly  and  in  the  integral  on  the  right  side 
of  Eq.  (22)  can  be  replaced  with  the  surface  temperature  and  the  tempera- 
ture from  the  temperature  profile  in  the  atmosphere,  respectively.  For 
microwave  remote  sensing,  much  of  the  lower  boundary  (water)  has  an  emis- 
sivlty  e significantly  less  than  unity.  Consequently,  two  changes  must 
be  Introduced  in  the  radiance  expression:  the  surface  emission  must  in- 
clude the  emissivity,  and  a reflection  (1  - eg)  of  downwelling  radiation 
at  the  surface  must  be  introduced.  Thus  the  expression  for  microwave 
measurements  in  Kelvins  at  a given  frequency  becomes  (cf.  Grody)17 

' dl  (,)  fS  dT  (2) 

T<*>  -3T-  + (l  - ‘s>  Vs)  L ™ -i — 

J H 

(24) 

where  t„(z)  is  always  measured  along  a path  extending  from  H to  z,  t_(z) 
n o 

is  always  referenced  at  the  lower  boundary  S,  and  th(S)  ■ Tg(H).  Strictly 
speaking,  with  the  introduction  of  an  emissivity  less  than  one,  the  true 


*sTs 


vs> 


/, 


52 


► 


I 


i 


weighting  functions  are  no  longer  given  by  dT^(z)/dz.  Instead,  the 
weighting  functions  are  now  dependent  on  the  emissivity  and  on  the  prod- 
uct T„(S)  dT0(z)/dz.  This  complication  was  not  introduced  in  the  previous 

H S 

illustrations. 

B.  Infrared  Data 

For  the  infrared  computations,  the  lower  boundary  was  assumed  to  be 
a blackbody.  Computations  were  made  for  all  20  cm  1 spectral  intervals 
involved  in  the  matching  of  weighting  functions.  Both  the  Midlatitude 
Summer  and  Midlatitude  Winter  soundings  were  applied,  and  compuations  at 
0°  zenith  and  40°  zenith  angles  were  included  for  both  soundings.  In 
addition,  without  altering  the  soundings,  computations  were  performed  with 
the  blackbody  lower  boundary  elevated  to  3,  6,  and  9 km.  To  facilitate 
comparisons,  each  computed  radiance  was  converted  to  an  equivalent  Planck- 
ian  temperature  through  the  inverse  to  Eq.  (23)  with  the  computed  radiance 
(R)  replacing  the  Planck  radiance: 

Tb(z)  * C2v1  j loge[(C1vi3/R)  + 1]  . (25) 

The  height  profiles  of  pressure  and  temperature  defined  for  the  two 
AFGL  soundings  are  listed  up  to  100  km  in  Table  5 (densities  of  air, 
water  vapor,  and  ozone  are  not  repeated  here).  These  data  can  be  compared 
with  the  computed  radiative  temperatures.  Tables  6a  and  6b  list  the  com- 
puted radiances  and  equivalent  temperatures,  grouped  according  to  antici- 
pated matches  in  the  weighting  functions.  These  data  constitute  a good 
set  of  background  information  (molecular  atmosphere)  for  evaluation  of 
relative  sensing  and  for  extended  study  with  the  introduction  of  cloud 
information.  Some  Initial  assessments  can  be  obtained  from  simple  mani- 
pulations of  these  data,  as  mentioned  in  the  subsequent  discussion  of 
Subsection  D. 

| 

Inspection  of  the  radiance  tables  reveals  that  the  excellent  match 
between  445-465  cm  and  1365-1385  cm  1 is  maintained  for  both  viewing 
angles  in  both  seasons.  Water  vapor  dominates  both  spectral  channels, 
as  it  does  for  the  intervals  530-550  cm'1  and  1925-1945  cm  1.  As  antl- 
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Table  5 


HEIGHT  PROFILES  OF  PRESSURE  AND  TEMPERATURE  FOR 
MIDLATITUDE  SUMMER  AND  MIDLATITUDE  HUNTER 


Height 

(km) 

Summer 

Winter 

Pressure 

(mb ) 

Temperature 

00 

Pressure 

(mb) 

Temperature 

OP 

0 

1013 

294 

1018.0 

272.2 

1 

902 

290 

897.3 

268.7 

2 

802 

285 

789.7 

265.2 

3 

710 

279 

693.8 

261.7 

4 

628 

273 

608.1 

255.7 

5 

554 

267 

531.3 

249.7 

6 

487 

261 

462.7 

243.7 

7 

426 

255 

401.6 

237.7 

8 

372 

248 

347.3 

231.7 

9 

324 

242 

299.2 

225.7 

10 

281 

235 

256.8 

219.7 

11 

243 

229 

219.9 

219.2 

12 

209 

222 

188.2 

219.7 

13 

179 

216 

161.0 

218.2 

14 

153 

216 

137.3 

217.7 

15 

130.0 

216 

117.8 

217.2 

16 

111.0 

216 

100.7 

216.7 

17 

95.0 

216 

86.1 

216.2 

18 

81.2 

216 

73.5 

215.7 

19 

69.5 

217 

62.8 

215.2 

20 

69.5 

218 

53.7 

216.2 

21 

51.0 

219 

45.8 

215.2 

22 

43.7 

220 

39.1 

215.2 

23 

37.6 

222 

33.4 

215.2 

24 

32.2 

223 

28.6 

215.2 

25 

27.7 

224 

24.3 

215.2 

30 

13.2 

234 

11.1 

217.4 

36 

6.52 

246 

5.18 

227.8 

40 

3.33 

258 

2.53 

243.2 

45 

1.76 

270 

1.29 

258.6 

60 

.9510 

276 

.6820 

265.7 

70 

.0671 

213 

.0467 

230.7 

100 

.0003 

210 

.0003 

210.2 

Table  6 


COMPUTED  INFRARED  RADIANCES  AND  EQUIVALENT  TEMPERATURES 
FOR  TWO  ANGLES  OVER  FOUR  LOWER  BOUNDARY  HEIGHTS 
(a)  Midlatitude  Summer 


Spectral  Surfe.e 

0-Degree 

Zenith  Angle 

40-Degree 

Zenith  Angl< 

Interval  Height 
(cm-1)  (kn>) 

Radiance 
mW/ (m^sr/ cm 

Temperature 
) (K) 

Radiance 
mW/ (m^sr/cm) 

Temperature 

00 

700-720 

0 

67.450 

245.45 

64.055 

242.49 

3 

66.484 

244.62 

63.363 

241.87 

6 

63.219 

241.74 

60.956 

239.71 

9 

55.309 

234.43 

54.368 

233.52 

2240-2260 

0 

0.286 

247.70 

0.232 

243.82 

3 

0.273 

246.80 

0.226 

243.29 

6 

0.223 

243.10 

0.195 

240.67 

9 

0.138 

233.59 

0.130 

233.57 

720-740 

0 

85.527 

261.90 

81.763 

259.04 

3 

80.529 

258.09 

77.789 

255.95 

6 

69.219 

248.98 

67.664 

247.67 

9 

55.645 

236.88 

55.058 

236.33 

725-745 

0 

84.302 

261.46 

79.883 

258.08 

3 

80.921 

258.88 

77.364 

256.11 

6 

70.971 

250.95 

69.131 

249.42 

9 

56.478 

238.21 

55.858 

237.63 

735-755 

0 

90.840 

267.23 

86.627 

264.16 

3 

85.436 

263.28 

82.385 

260.98 

6 

71.895 

252.73 

70.339 

251.45 

9 

55.662 

238.52 

55.087 

237.98 

2180-2200 

0 

1.545 

278.81 

1.395 

276.31 

3 

1.070 

270.03 

0.999 

268.46 

6 

0.571 

256.22 

0.548 

255.39 

9 

0.247 

239.92 

0.243 

239.57 

880-900 

0 

106.179 

292.16 

105.377 

291.67 

3 

85.894 

278.83 

86.820 

278.78 

6 

62.577 

260.99 

62.562 

260.98 

9 

42.480 

242.00 

42.480 

242.00 

1190-1210 

0 

57.018 

293.06 

56.709 

292.79 

3 

42.260 

278.92 

42.236 

278.89 

6 

27.605 

260.99 

27.604 

260.99 

9 

16.416 

242.00 

16.416 

242.00 

2500-2520 

0 

0.829 

292.81 

0.816 

292.45 

3 

0.439 

278.48 

0.436 

278.32 

6 

0.182 

260.78 

0.182 

260.72 

9 

0.062 

241.92 

0.062 

241.89 
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Table  6(a)  (continued) 


Spectral  Surface  0-Degree  Zenith  Anglr  40-Pegree  Zenith  Angle 

Interval  Height  Radiance  Temperature  Radiance  Temperature 

(cm~*)  (km)  mW/ (m^sr/cm)  (K) mW/ (m^sr/cm)  (K)  . 


445-465 

0 

100.248 

261.79 

98.325 

259.94 

3 

97.356 

259.00 

95.569 

257.27 

6 

94.415 

256.14 

93.539 

255.28 

9 

79.780 

241.38 

79.687 

241.28 

1365-1385 

0 

16.162 

261.75 

15.242 

259.74 

3 

15.520 

260.35 

14.703 

258.51 

6 

13.375 

255.36 

13.023 

254.48 

9 

8.416 

240.95 

8.372 

240.80 

530-550 

0 

115.700 

273.05 

11-.323 

271.19 

3 

113.313 

271.18 

111.058 

269.40 

6 

98.531 

259.21 

98.199 

258.93 

9 

78.499 

241.70 

78.454 

241.66 

1925-1945 

0 

3.858 

277.98 

3.620 

276.23 

3 

3.185 

272.76 

3.072 

271.80 

6 

1.868 

259.21 

1.847 

258.95 

9 

0.855 

241.65 

0.853 

241.59 

530-540 

0 

119.770 

276.23 

117.449 

274.42 

3 

117.220 

274.24 

115.392 

272.80 

6 

99.472 

259.89 

99.274 

259.72 

9 

78.865 

241.82 

78.838 

241.79 

535-545 

0 

116.862 

273.96 

114.467 

272.09 

3 

114.466 

272.09 

112.297 

270.38 

6 

98.715 

259.36 

98.416 

259.11 

9 

78.510 

241.71 

78.466 

241.67 

540-550 

0 

111.935 

270.21 

109.424 

268.24 

3 

109.556 

268.34 

106.953 

226.28 

6 

97.222 

258.40 

96.727 

257.99 

9 

77.785 

241.55 

77.715 

241.48 

1925-1935 

0 

4.524 

281.97 

4.279 

280.39 

3 

3.518 

274.95 

3.433 

274.28 

6 

1.955 

259.83 

1.941 

259.66 

9 

0.880 

241.77 

0.878 

241.74 

1930-1940 

0 

4.387 

281.60 

4.144 

279.98 

3 

3.435 

274.80 

3.348 

274.11 

6 

1.913 

259.80 

1.899 

259.62 

9 

0.860 

241.76 

0.859 

241.73 

1935-1945 

0 

3.258 

273.87 

3.039 

272.01 

3 

2.818 

270.02 

2.687 

268.79 

6 

1.763 

258.21 

1.735 

157.92 

9 

0.830 

241.46 

0.827 

241.39 
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Table  6 (continued) 
(b)  Midlatitude  Winter 


SgfigUftl  Surface 

fl-ngarg«  Zonlth  Angle 

40-DeRree 

Zenith  Angle 

Interval  Height 
(c”"1)  (km) 

Radiance 

Temperature 

(K) 

Radiance 
mW/(ni  sr/cm) 

Temperature 

(K) 

700-720 

0 

57.579 

236.60 

54.890 

234.03 

3 

56.482 

235.55 

54.049 

233.21 

6 

52.138 

231.33 

50.712 

229.91 

9 

44.040 

222.91 

43.714 

222.55 

2240-2260 

0 

0.156 

236.70 

0.129 

233.45 

3 

0.148 

235.79 

0.124 

232.85 

6 

0.112 

231.07 

0.101 

229.37 

9 

0.066 

222.61 

0.064 

222.21 

720-740 

0 

69.597 

249.30 

66.988 

247.09 

3 

65.005 

245.38 

63.043 

243.66 

6 

54.242 

235.54 

53.282 

234.62 

9 

42.704 

223.67 

42.490 

223.43 

725-745 

0 

69.742 

249.93 

66.529 

247.21 

3 

66.262 

246.98 

63.788 

244.84 

6 

55.533 

237.32 

54.478 

236.31 

9 

42.739 

224.27 

42.533 

224.05 

735-755 

0 

73.879 

254.34 

70.968 

251.97 

3 

68.481 

249.91 

66.390 

248.14 

6 

55.259 

238.15 

54.345 

237.27 

9 

41.752 

224.31 

41.550 

224.08 

2180-2200 

0 

0.722 

261.23 

0.663 

259.38 

3 

0.527 

254.58 

0.496 

253.33 

6 

0.252 

240.26 

0.244 

239.67 

9 

0.102 

224.79 

0.101 

224.64 

880-900 

0 

76.378 

271.94 

76.275 

271.86 

3 

63.382 

261.66 

63.369 

261.65 

6 

44.085 

243.70 

44.085 

243.70 

9 

28.942 

225.70 

28.942 

225.70 

1190-1210 

0 

36.156 

272.07 

36.124 

272.04 

3 

28.091 

261.68 

28.087 

261.68 

6 

17.254 

243.70 

17.254 

243.70 

9 

9.803 

225.70 

9.803 

225.70 

2500-2520 

0 

0.313 

271.36 

0.309 

271.10 

3 

0.187 

261.27 

0.186 

261.14 

6 

0.068 

243.54 

0.068 

243.49 

9 

0.021 

225.66 

0.021 

225.65 

i 

i 
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Table  6 (b)  (concluded) 


Spectral  Surface 

0-Degree 

Zenith  Angle 

40-Degree  Zenith  Angle 

Interval  Height 
(cm-^)  (km) 

Radiance 
mW/ (m^sr/cm 

Temperature 
) . 00 

Radiance 
mW/ (m^sr/ cm) 

Temperature 

00 

445-465 

0 

93.209 

254.96 

91.694 

253.46 

3 

91.186 

252.96 

89.696 

251.48 

6 

81.055 

242.70 

80.885 

242.52 

9 

65.259 

225.67 

65.255 

225.67 

1365-1385 

0 

12.987 

254.39 

12.364 

252.79 

3 

12.250 

252.49 

11.746 

251.15 

6 

8.704 

241.94 

8.625 

241.68 

9 

4.792 

225.49 

4.786 

225.46 

530-550 

0 

103.188 

263.06 

101.721 

261.85 

3 

97.906 

258.69 

97.260 

258.15 

6 

80.191 

243.25 

80.123 

243.19 

9 

61.938 

225.67 

61.934 

225.66 

1925-1945 

0 

2.367 

265.06 

2.279 

264.11 

3 

1.838 

258.82 

1.805 

258.39 

6 

0.917 

243.11 

0.913 

243.03 

9 

0.377 

225.62 

0.377 

225.61 

530-540 

0 

105.932 

265.22 

104.645 

264.17 

3 

99.461 

259.88 

99.067 

259.55 

6 

80.624 

243.43 

80.588 

243.40 

9 

62.263 

225.68 

62.261 

225.68 

535-545 

0 

103.871 

263.61 

102.430 

262.43 

3 

98.261 

258.98 

97.678 

258.50 

6 

80.208 

243.27 

80.145 

243.21 

9 

61.935 

225.66 

61.930 

225.66 

540-550 

0 

100.200 

260.84 

98.527 

259.47 

3 

95.927 

257.32 

94.992 

256.54 

6 

79.382 

243.00 

79.270 

242.89 

9 

61.253 

225.64 

61.245 

225.64 

1925-1935 

0 

2.627 

267.22 

2.552 

266.48 

3 

1.954 

259.83 

1.932 

259.54 

6 

0.946 

243.32 

0.944 

243.26 

9 

0.387 

225.65 

0.387 

225.64 

1930-1940 

0 

2.561 

267.06 

2.485 

266.29 

3 

1.911 

259.77 

1.888 

259.47 

6 

0.926 

243.31 

0.923 

243.25 

9 

0.378 

225.65 

0.378 

225.64 

1935-1945 

0 

2.097 

262.53 

2.003 

261.39 

3 

1.699 

257.43 

1.656 

256.81 

6 

0.884 

242.81 

0.879 

242.68 

9 

0.368 

225.58 

0.367 

225.56 
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cipated  from  the  weighting  function  profiles  discussed  in  Section  VI,  the 
540  cm  ^ spectral  region  is  associated  with  stronger  attenuation  and, 
therefore,  lower  temperatures  than  the  1935  cm  ^ region.  The  discrepancy 
in  temperature  is  reduced  in  the  winter  because  of  the  reduced  water 
vapor.  In  the  summer  the  discrepancy  is  maintained  above  3 km  but  in 
* winter  the  discrepancy  is  confined  below  3 km.  Much  of  the  difference 

between  channels  is  due  to  the  much  stronger  continuum  at  540  cm 


t 


A review  of  the  10  cm  intervals  shows  that  for  both  the  540  and 
1935  cm  ^ regions  the  absorption  (and  weighting  functions)  showed  an  in- 
crease with  wavenumber.  The  computed  radiances  show  that  the  first  two 
subdivisions  of  the  1935  cm  ^ interval  are  nearly  the  same,  with  a rapid 
reduction  for  the  last  subdivision.  In  contrast,  the  540  cm  subdivi- 
sions change  more  uniformly,  and  535-545  is  fairly  representative  of  the 
larger  530-550  cm  ^ interval.  (Of  course,  the  average  of  the  intervals 
530-540  and  540-550  is  closest  to  the  530-550  result.)  In  this  case 
changing  the  size  of  the  interval  has  little  impact  on  the  rad'ance,  but 
the  same  is  not  true  for  the  1935  cm  ^ interval,  where  shrinking  the  in- 
terval increases  the  radiance. 


The  infrared  results  discussed  thus  far  have  been  based  on  the  trans- 
mittance model  relating  the  transmittance  function  i|i  to  A/VP . With  the 

exception  of  continuum  transmittance,  models  for  other  spectral  intervals 
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relate  V to  AW/[(W/P)  + B ] ‘ . From  Tables  6a  and  6b  it  is  apparent 
that  the  700-720  cm  * interval  is  well  matched  with  the  2240-2260  cm 
interval,  although  equivalent  temperatures  are  a bit  lower  for  the  summer 
zenith  position.  On  the  other  hand,  the  anticipated  match  between  720- 
740  cm  * and  2180-2200  cm  ^ is  not  there;  the  720-740  cm  interval  is 
associated  with  greater  attenuation  in  the  upper  atmosphere  and  lower 
equivalent  temperatures.  Furthermore,  the  725-745  cm  interval  (which 
yields  results  very  similar  to  the  720-740  cm  * interval)  and  the  735-755 
cm  Interval  also  remain  colder  than  the  2180-2200  cm  * interval.  To 
find  a good  match  requires  a more  careful  surveillance  of  both  spectral 
regions. 


In  each  of  the  three  window  intervals  centered  at  890,  1200,  and 
2510  cm  ^ , only  the  continuum  was  considered.  The  1200  cm  interval 
demonstrated  the  greatest  attenuation,  but  in  winter  the  decrease  in 
water  vapor  meant  that  the  2510  cm  ^ interval  with  the  nitrogen  continuum 
retained  slightly  higher  attenuation  than  the  890  cm*1'  interval.  Empiri- 
cal checks  are  needed  for  these  channels,  since  it  appears  that  all  of 
the  window  channels  are  too  transparent;  greater  attenuation  was  antici- 
pated. 

C.  Microwave  Data 

All  of  the  monochromatic  microwave  computations  were  performed  for 
the  same  two  angles  for  each  of  the  same  t\  ? summer  and  winter)  sound- 
ings, with  the  lower  boundary  taken  successively  at  0,  3,  6,  and  9 km. 
However,  two  additional  types  of  computations  were  made  for  each  frequency. 
The  emissivity  of  the  lower  boundary  was  set  at  1.0,  0.5,  and  0.25.  Next, 
two  types  of  surface  reflectance  were  considered:  specular  and  diffuse. 

For  the  specular  reflectance,  the  downwelling  radiance  at  the  lower  boun- 
dary was  computed  only  at  the  same  zenith  angle  as  the  emergent  beam 
(cf.  Grody).17  For  the  diffuse  reflectance,  Instead  of  multiplying  the 
vertical  optical  path  by  the  secant  of  the  zenith  angle,  a factor  of  1.6 
was  used  regardless  of  zenith  angle.  This  factor  approximates  the  trans- 
formation from  downward  beam  to  flux  transmittance  (the  main  concern  here 
was  the  magnitude  of  the  departure  from  the  specular  case).  A major  dif- 
ference is  that  for  a comparison  of  radiances  along  the  two  slant  paths, 
the  downwelling  radiance  for  the  diffuse  case  is  the  same  for  both  0° 
and  40°  zenith  angles.  This  should  reduce  any  limb  brightening,  but 
could  enhance  limb  darkening. 

Tables  7a  and  7b  present  results  for  four  frequencies  in  the  oxygen 
band,  for  summer  and  winter;  Tables  7c  and  7d  cover  the  results  at  four 
frequencies  near  the  183  GHz  line  of  water  vapor,  and  Tables  7e  and  7f 
cover  several  window  frequencies  as  well  as  the  weak  ’’2.235  GHz  line  of 
water  vapor.  Presumably,  the  emissivity  of  1.0  would  give  results  most 
consistent  with  viewing  over  land  while  the  < * 0.5  might  be  valid  for 
the  sea  surface  on  the  average. 


60 


Table  7 

COMPUTED  MICROWAVE  BRIGHTNESS  TEMPERATURES  FOR  TWO  ZENITH  ANGLES 
AND  VARIED  SURFACE  EUISSIVITY,  HEIGHT,  AND  REFLECTANCE 


(a)  Oxygen,  Midlatitude  Summer 


Lmwr  cuDegree  Zenith  Angle  aruPegree  Zenith  Angle 

Boundary  Surfaoe  Qaisairlty  Surface  SmissiTity 


(k») 

1*0 

0*6 

0.26 

1*0 

0.6 

0.26 

Specular 

0 

271.90 

267.69 

52 

250.43 

.85  GHz 

267.31 

260.08 

256.46 

3 

266.69 

226.84 

205.91 

262.52 

235.02 

221.28 

6 

254.64 

185.91 

151.59 

252.84 

196.18 

167.85 

9 

239.6*' 

150.98 

106.68 

238.90 

168.19 

117.84 

Diffuse 

0 

271.90 

264.16 

260.29 

267.31 

261.96 

259.29 

3 

265.69 

237.44 

223.32 

262.52 

239.34 

227.76 

6 

254.54 

197.28 

168.65 

252.84 

201.07 

175.19 

9 

239.57 

158.57 

118.07 

238.90 

161.72 

123.14 

Specular 

0 

261.01 

257.29 

_53 

255.43 

.5. 

GHz 

255.40 

254.09 

253.44 

3 

267.46 

240.25 

231.65 

253.07 

243.73 

239.06 

6 

249.57 

207.99 

187.20 

246.88 

217.25 

202.44 

9 

237.36 

171.21 

138.14 

236.15 

180.92 

153.30 

Diffuse 

0 

261.01 

259.44 

258.65 

255.40 

254.53 

254.10 

3 

267.46 

247.46 

242.47 

253.07 

245.91 

242.34 

6 

249.57 

219.19 

204.00 

246.88 

221.52 

208.84 

9 

237.35 

181.63 

163.77 

236.15 

185.44 

160.08 

Specular 

0 

261.11 

250.36 

54 

249.98 

• 0. 

GHz 

245.31 

245.13 

245.04 

3 

249.35 

245.36 

240.37 

244.36 

241.92 

240.70 

6 

244.30 

222.56 

211.69 

240.88 

228.02 

221.59 

9 

234.88 

189.22 

166.39 

233.22 

199.03 

181.94 

Diffuse 

0 

261.11 

260.83 

260.68 

246.31 

245.19 

245.13 

3 

249.35 

246.49 

245.06 

244.36 

242.66 

241.79 

6 

244.30 

230.74 

223.96 

240.88 

230.69 

225.59 

9 

234.88 

200.10 

182.70 

233.22 

203.35 

188.42 

Specular 

0 

234.01 

233.99 

54 

238.99 

.9_ 

GHz 

229.66 

229.66 

229.56 

3 

233.81 

238.56 

233.43 

229.50 

229.46 

229.42 

6 

232.64 

230.04 

228.74 

226.96 

228.08 

227.64 

9 

228.72 

215.64 

209.10 

226.46 

219.62 

216.19 

Diffuse 

0 

234.01 

234.00 

234.00 

229.56 

229.56 

229.66 

3 

233.81 

233.70 

233.64 

229.60 

229.46 

229.44 

6 

232.64 

231.49 

230.91 

228.96 

228.36 

228.05 

9 

228.72 

221.36 

217.66 

226.48 

221.30 

218.70 
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Table  7 (continued) 

(b)  Oxygen,  Midlatitude  Winter 

Uver  Q-Pegree  Zenith  Angle  40-Degree  Zenith  Angle 


Bourn  try  Surface  Bbiiaaielty  Surface  fimissiTity 


(km/ 

1.0 

0.5 

0.25 

1.0 

0.5 

0.25 

Specular 

0 

255.13 

242.88 

52_ 

236.76 

.85  GHz 

251.58 

245.54 

242.62 

3 

250.86 

211.99 

192.56 

248.29 

221.18 

207.62 

6 

239.39 

172.11 

138.47 

238.26 

181.91 

163.74 

9 

224.85 

139.09 

96.21 

224.61 

145.63 

106.13 

Diffuse 

0 

255.13 

248.72 

245.52 

251.58 

247.19 

245.00 

3 

250.86 

223.08 

209.19 

248.29 

225.34 

213.87 

6 

239.39 

182.61 

154.22 

238.26 

186.49 

160.60 

9 

224.85 

146.76 

106.20 

224.61 

148.75 

110.82 

Specular 

0 

247.13 

243.86 

53. 

242.22 

;5_ 

GHz 

242.83 

241.71 

241.14 

3 

244.62 

226.74 

217.80 

241.18 

231.19 

226.20 

6 

236.23 

193.25 

171.77 

234.46 

202.93 

187.17 

9 

224.04 

157.85 

124.75 

223.61 

167.18 

138.97 

Diffuse 

0 

247.13 

245.81 

245.15 

242.83 

242.10 

241.74 

3 

244.62 

234.03 

228.74 

241.18 

233.46 

229.60 

6 

236.23 

204.12 

188.07 

234.46 

207.19 

193.55 

9 

224.04 

167.42 

139.12 

223.61 

171.42 

145.32 

Specular 

0 

239.70 

239.04 

54 

238.70 

.0_ 

GHz 

235.29 

235.13 

235.06 

3 

238.43 

231.83 

226.53 

234.60 

231.79 

230.38 

e 

232.84 

209.06 

197.16 

230.63 

215.92 

208.57 

9 

223.17 

176.64 

151.88 

222.56 

185.86 

167.50 

Diffuse 

0 

239.70 

239.46 

239.35 

235.29 

236.19 

235.14 

3 

238.43 

235.21 

233.60 

234.60 

232.60 

231.61 

6 

232.84 

217.57 

209.93 

230.63 

218.81 

212.90 

9 

223.17 

186.21 

167.74 

222.56 

190.18 

173.99 

Specular 

0 

226.93 

226.93 

54 

226.92 

.9_ 

GHz 

223.81 

223.81 

223.81 

3 

226.79 

226.49 

226.34 

223.77 

223.71 

223.67 

6 

225.38 

222.18 

220.58 

223.09 

221.93 

221.35 

9 

220.99 

206.76 

198.14 

220.16 

211.73 

207.52 

Diffuse 

0 

226.93 

226.93 

226.93 

223.81 

223.81 

223.81 

3 

226.79 

226.66 

226.69 

223.77 

223.72 

223.70 

6 

225.38 

223.95 

223.24 

223.09 

222.29 

221.90 

9 

220.99 

212.30 

207.98 

220.16 

213.78 

210.58 
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Table  T (continued) 


4 

(c)  Water  Vapor,  Midlatitude  Summer 

Lower  0-Degree  Zenith  Angle  40-Degree  Zenith  Angle 

* Boundary  Surface  BmisslTity  Surface  Bnissivlty 


(km) 

1.0 

0.5 

0.25 

1.0 

0.5 

0.25 

Specular 

0 

286.59 

271.76 

162. 

264.34 

5 GHz 
284.98 

277.48 

273.72 

3 

277.30 

180.55 

132.18 

276.82 

190.12 

146.77 

6 

260.75 

137.13 

75.32 

260.67 

139.05 

78.23 

9 

241.98 

121.74 

61.62 

241.97 

121.96 

61.96 

Diffuse 

0 

286.59 

278.89 

275.04 

284.98 

279.56 

276.85 

3 

277.30 

190.41 

146.96 

276.82 

194.65 

153.57 

6 

260.76 

139.08 

78.25 

260.67 

140.01 

79.68 

9 

241.98 

121.97 

61.96 

241.97 

122.07 

62.13 

Specular 

0 

283.58 

279.68 

170 

277.72 

GHz 

281.63 

280.28 

279.60 

3 

276.36 

198.58 

159.69 

275.64 

210.34 

177.70 

6 

260.60 

140.97 

81.16 

260.48 

143.93 

85.66 

9 

241.96 

122.19 

62.31 

241.95 

122.55 

62.85 

Diffuse 

0 

283.58 

282.12 

281.39 

281.63 

280.79 

280.38 

3 

276.36 

210.81 

178.03 

275.64 

215.69 

185.71 

6 

260.60 

143.99 

85.69 

260.48 

145.41 

87.88 

9 

241.96 

122.55 

62.85 

241.95 

122.73 

63.12 

Specular 

0 

272.17 

272.17 

178 

272.17 

GHz 

269.71 

269.71 

269.71 

3 

270.47 

255.49 

247.99 

268.66 

260.92 

267.05 

6 

259.41 

167.52 

121.68 

268.97 

176.24 

134.88 

9 

241.86 

125.87 

67.88 

241.80 

127.29 

70.03 

Diffuse 

0 

272.17 

272.17 

272.17 

269.71 

269.71 

269.71 

3 

270.47 

262.48 

268.48 

268.66 

262.99 

260.16 

6 

259.41 

176.61 

135.06 

258.97 

180.39 

141.10 

9 

241.85 

127.31 

70.04 

241.80 

128.00 

71.10 

Specular 

0 

262.13 

262.13 

182 

262.13 

GHz 

249.53 

249.55 

249.53 

3 

262.13 

252.13 

252.13 

249.63 

249.53 

249.53 

6 

260.69 

242.26 

238.03 

248.73 

244.95 

243.06 

9 

240.47 

161.53 

122.06 

240.05 

170.69 

136.86 

Diffuse 

0 

252.13 

252.13 

252.13 

249.53 

249.53 

249.53 

3 

252.13 

262.13 

252.13 

249.53 

249.53 

249.63 

6 

260.69 

246.69 

244.69 

248.73 

246.11 

244.80 

9 

240.47 

170.84 

136.02 

240.05 

174.81 

142.19 
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Table  T (continued) 


(d)  Water  Vapor,  Midlatitude  Winter 

Lower  0-Degree  Zenith  Angle  40~Pegree  Zenith  Angle 

Boundary  Surfaoe  Eroissirity  Surface  I&iissiTity 


(km) 

1.0 

0.5 

0.25 

1.0 

0.5 

0.25 

Specular 

0 

269.76 

206.23 

162.5 

174.46 

GHz 

269.12 

218.54 

193.25 

3 

261.11 

147.34 

90.46 

260.94 

151.84 

97.29 

6 

243.66 

223.66 

63.51 

243.63 

124.07 

64.26 

9 

225.70 

112.99 

56.64 

225.70 

113.03 

66.70 

Diffuse 

0 

269.76 

218.98 

193.59 

269.12 

223.89 

201.28 

3 

261.11 

151.93 

97.35 

260.94 

154.07 

100.63 

6 

243.65 

124.07 

64.29 

243.63 

124.32 

64.67 

9 

225.70 

113.03 

56.70 

225.70 

113.05 

66.73 

Specular 

0 

268.61 

227.79 

170 

207.43 

GHz 

267.60 

239.21 

225.02 

3 

260.76 

156.14 

103.82 

260.50 

162.65 

113.73 

6 

243.61 

124.59 

65.07 

243.59 

125.39 

66.30 

9 

225.69 

113.08 

56.77 

225.69 

113.15 

56.87 

Diffuse 

0 

268.61 

239.90 

225.60 

267.60 

243.80 

231.89 

3 

260.76 

162.80 

113.82 

260.50 

165.82 

118.48 

6 

243.61 

125.41 

66.30 

243.59 

125.80 

66.91 

9 

225.69 

113.15 

56.87 

225.69 

113.18 

56.93 

Specular 

0 

262.14 

260.16 

178 

259.17 

GHz 

260.40 

259.82 

259.53 

3 

258.34 

202.79 

175.02 

257.48 

214.48 

192.99 

6 

243.36 

132.51 

77.09 

243.25 

135.51 

81.64 

9 

225.68 

113.80 

67.87 

225.67 

114.09 

58.29 

Diffuse 

0 

262.14 

261.48 

261.16 

260.40 

260.06 

259.89 

3 

268.34 

215.08 

193.45 

257.48 

219.63 

200.55 

6 

243.36 

136.67 

81.87 

243.25 

137.01 

83.88 

9 

226.68 

114.09 

58.30 

226.67 

114.23 

58.51 

Specular 

0 

247.77 

247.77 

182 

247.77 

GHz 

246.46 

245.46 

245.45 

3 

247.53 

246.24 

245.59 

245.36 

245.00 

244.82 

6 

240.89 

182.67 

153.66 

240.16 

193.48 

170.13 

9 

226.43 

122.94 

71.70 

226.36 

126.78 

76.00 

Diffuse 

0 

247.77 

247.77 

247.77 

245.46 

245.46 

245.45 

3 

247.63 

247.09 

246.86 

245.36 

245.14 

245.03 

6 

240.89 

193.97 

170.61 

240.16 

198.24 

177.28 

9 

225.43 

126.82 

76.02 

225.36 

127.19 

78.11 
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Table  f (oontinued) 


r 

1 (e)  Cxygen-Water  Vapor,  Midlatitude  Summer 


Lower 

0-Degree  Zenith  Angle 

40-Degre 

:e  Zenith  Angle 

• 

Boundary 

Surfaoe  BnlssWity 

Surface  Qnis« 

ilrity 

(km) 

1.0 

0.6 

0,25 

1.0 

0.5 

0.25 

# 

Specular 

19.35  GHz 

0 

290.16 

190.88 

141.23 

289.08 

200.80 

156.66 

5 

277.02 

159.23 

100.34 

276.45 

164.48 

108.49 

6 

260.11 

139.79 

79.63 

259.85 

142.42 

83.70 

9 

241.68 

125.35 

67.19 

241.58 

126.63 

69.15 

Diffuse 

0 

290.16 

201.47 

157.13 

289.08 

205.64 

163.91 

3 

277.02 

164.79 

106.66 

276.45 

167.15 

112.51 

S 

260.11 

142.56 

83.78 

259.85 

143.78 

85.74 

9 

241.68 

126.67 

69.17 

241.58 

127.28 

70.13 

Specular 

22.235  GHz 

0 

288.63 

210.85 

171.97 

287.17 

222.88 

190.73 

3 

276.43 

167.92 

113.66 

275.69 

175.04 

124.72 

6 

259.94 

142.77 

84.19 

259.62 

146.18 

89.46 

9 

241.66 

126.39 

68.76 

241.55 

127.96 

T1.16 

Diffuse 

0 

288.63 

223.87 

191.49 

287.17 

228.49 

199.15 

3 

276.43 

175.47 

124.99 

275.69 

178.61 

130.06 

6 

259.94 

146.35 

89.55 

259.62 

147.93 

92.06 

9 

241.66 

128.01 

71.18 

241.55 

128.76 

72.36 

Specular 

31.4 

GHz 

0 

290.01 

190.93 

141.40 

288.88 

200.84 

156.83 

3 

276.90 

159.86 

101.34 

276.29 

165.24 

109.72 

6 

260.04 

140.38 

80.55 

259.76 

143.17 

84.87 

9 

241.65 

125.70 

67.72 

241.55 

127.07 

69.83 

Diffuse 

0 

290.01 

201.55 

157.32 

288.88 

205.69 

164.09 

3 

276.90 

165.58 

109.92 

276.29 

167.99 

113.84 

6 

260.04 

143.52 

84.95 

259.76 

144.61 

87.03 

9 

241.65 

127.12 

69.85 

241.55 

127.77 

70.88 

Specular 

37.0 

GHz 

0 

269.67 

193.61 

145.58 

288.46 

203.88 

181.59 

. 

3 

276.71 

161.49 

103.88 

276.06 

167.24 

112.84 

6 

259.96 

141.25 

81.90 

269.66 

144.27 

86.58 

9 

241.62 

126.12 

68.38 

241.60 

127.61 

70.67 

Diffuse 

0 

299.67 

204.66 

162.14 

288.46 

206.88 

169.09 

3 

276. 7i 

167.61 

113.06 

276.06 

170.17 

117.23 

1 

6 

259.96 

144.43 

86.67 

259.66 

145.82 

88.91 

9 

241.62 

127.67 

70.70 

241.60 

128.38 

71.81 

l 
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Table  7 (oonoluded) 


(f)  Oxygen-Water  Vapor,  Midlatitude  Winter 


Lower 

0-Degree 

Zenith 

Angle 

40-Degree  Zenith  Angle 

Boundary 

Surface  Beissi.Ti.ty 

Surfaoe  Emissirity 

0-) 

1.0 

0.6 

0.26 

1.0 

0.5 

0.25 

Specular 

19.35  GHz 

0 

269.28 

172.71 

124.42 

268.45 

181.31 

137.74 

5 

260.02 

148.94 

93.40 

259.54 

153.77 

100.89 

6 

243.06 

130.46 

74.15 

242.90 

132.90 

77.90 

9 

226.58 

116.73 

62.30 

225.65 

117.86 

64.02 

Diffuse 

0 

269.28 

181.81 

138.07 

268.45 

185.51 

144.04 

3 

260.02 

154.04 

101.04 

259.54 

156.23 

104.57 

6 

243.08 

133.00 

77.96 

242.90 

134.15 

79.78 

9 

225.58 

117.88 

64.03 

225.55 

118.43 

64.88 

Specular 

22.235  GHz 

0 

268.79 

179.68 

135.13 

267.84 

189.34 

150.09 

3 

259.79 

152.19 

98.39 

259.24 

167.76 

107.02 

6 

243.02 

131.64 

75.95 

242.81 

134.40 

80.19 

9 

225.69 

117.48 

63.40 

225.56 

118.81 

65.43 

Diffuse 

0 

268.79 

189.94 

150.61 

267.84 

193.99 

167.06 

3 

259.79 

158.07 

107.20 

259.24 

160.57 

111.24 

6 

243.02 

134.50 

80.25 

242.81 

135.80 

82.30 

9 

225.59 

118.82 

65.44 

225.56 

119.48 

66.44 

Specular 

31.4 

GHz 

0 

269.06 

174.63 

127.26 

268.20 

183.42 

141.03 

3 

259.90 

150,07 

95.16 

259.37 

156.17 

103.06 

6 

243.03 

131.12 

75.17 

242.83 

133.76 

79.20 

9 

225.58 

117.06 

62.78 

225.54 

118.27 

64.64 

Diffuse 

0 

269.06 

183.96 

141.40 

268.20 

187.75 

147.63 

3 

259.90 

156.45 

103.23 

259.37 

167.75 

106.93 

6 

243.03 

133.85 

79.26 

242.83 

136.08 

81.20 

9 

225.58 

118.29 

64.65 

225.64 

118.89 

65.57 

Specular 

37.0 

i GHz 

0 

268.80 

177.34 

131.60 

267.86 

186.64 

146.04 

3 

269.73 

161.68 

97.66 

259.17 

167.13 

106.11 

6 

242.97 

131.95 

76.43 

242.763 

134.79 

80.80 

9 

225.56 

117.43 

63.36 

225.62 

118.76 

65.39 

Diffuse 

0 

268.80 

187.23 

146.46 

267.86 

191.16 

152.81 

3 

269.73 

167.44 

106.30 

259.17 

159.89 

110.26 

e 

242.97 

134.90 

80.86 

242.76 

136.23 

82.97 

9 

225.56 

118.78 

66.40 

225.62 

119.44 

66.39 
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The  results  for  the  oxygen  computations  show,  as  expected,  signifi- 
cant attenuation  that  increases  as  the  center  of  the  band  is  approached. 

The  lower  response  in  winter  reflects  the  lower  temperatures  for  that 
season.  As  required,  results  for  specular  and  diffuse  surface  reflectors 
are  identical  for  all  lower  boundary  altitudes  when  the  surface  emissivity 
is  1.0  (no  reflectance).  All  results  show  a limb  darkening  that  decreases 
as  the  boundary  altitude  increases.  With  the  surface  emissivity  set  at 
0.5,  and  looking  first  at  summer  results  with  the  lower  boundary  at  0 km, 
all  computations  show  a limb  darkening  except  for  the  52.85  GHz  frequency 
with  the  specular  surface.  Here  the  increased  downwelling  radiance  at 
40°  zenith  angle  (relative  to  the  0°  zenith  angle)  and  the  resulting  in- 
creased upward  reflected  radiance  with  less  attenuation  than  at  the  other 
frequencies,  produce  a limb  brightening.  It  also  is  apparent  that  the 
brightness  temperature  for  the  radiance  above  the  diffuse  surface  exceeds 
that  above  the  specular  surface,  especially  for  the  zenith  direction. 

Since  this  enhancement  results  from  the  larger  downwelling  radiance,  the 
effect  varies  with  the  altitude  of  the  surface  and  with  the  degree  of 
attenuation  for  the  given  frequency.  The  height  at  which  the  maximum  dif- 
ference occurs  increases  as  the  attenuation  increases.  A similar  sort  of 
pattern  prevails  for  the  surface  with  a 0.25  emissivity.  It  is  noteworthy 
that  with  the  surface  elevated  and  at  reduced  emissivity  the  emergent  radi- 
ances show  a limb  brightening  for  the  two  lowest  frequencies.  This  bright- 
ening is  most  pronounced  over  the  specular  surface  because  of  the  change 
m the  downwelling  radiance,  and  can  reach  substantial  magnitudes.  Pre- 
sumably, it  would  be  possible  to  test  empirically  which  type  of  reflecting 
surface  corresponds  most  closely  with  observation.  Because  most  elevated 
surfaces  would  not  be  likely  to  possess  an  emissivity  much  different  from 
unity,  comparisons  probably  should  be  based  only  on  surfaces  at  0 km  ele- 
vation. For  that  situation,  and  with  relatively  weak  atmospheric  atten- 
uation, significantly  reduced  limb  brightening  (or  enhanced  darkening) 
occurs  over  the  diffuse  surface. 

The  above  discussion  applies  also  to  the  winter  situation  and  is 
generally  consistent  with  the  results  obtained  as  the  183  GHz  line  of 
water  vapor  is  approached  (with  increasing  attenuation).  However,  the 
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brightness  remperatures  at  the  water  vapor  frequencies  generally  are 
higher  for  e = 1 and  for  all  e at  z = 0 in  summer,  with  a more  rapid  de- 
crease at  higher  boundary  elevations  when  e < 1.  In  winter,  with  less 
water  vapor,  temperatures  decrease  much  more  rapidly  than  for  0^  as  e is 
reduced,  as  well  as  when  the  height  of  the  boundary  is  increased. 

For  the  weak  22.235  water  vapor  line  and  three  nearby  window  regions, 
the  brightness  temperatures  are  summarized  for  the  same  conditions  as  be- 
fore, and  results  are  compatible  with  what  might  be  expected  for  weak  at- 
tenuation. When  the  surface  emissivity  is  reduced  to  0.5  or  0.25,  limb 
brightening  is  evident  everywhere  (even  with  the  0 km  boundary),  and  is 
smallest  over  the  diffuse  surface.  Differences  between  surface  types  are 
not  as  pronounced  as  with  some  of  the  previous  cases.  Although  window 
brightness  temperatures  are  high  for  the  blackbody  surface  at  0 km,  the 
temperature  decreases  more  rapidly  with  elevation  of  the  boundary  and 
with  reduction  of  the  emissivf-y  than  at  frequencies  with  greater  atten- 
uation. (High  boundaries  with  low  emissivities  are  not  realistic.) 

The  19.35  GHz  frequency  shows  the  strongest  return  over  a blackbody 
surface  (for  the  C^-H^O  combination)  but  results  are  not  much  different 
between  channels,  including  the  22.235  frequency.  For  the  surface  at  0-km 
altitude  the  brightness  temperature  over  the  diffuse  surface  is  distinctly 
larger. 

D.  Clouds 

* 

Without  computations  from  an  appropriate  physical  cloud  model  it  is 
not  possible  to  draw  inferences  concerning  cloud  physical  and  radiative 
properties  (cf.  Curran18)  from  the  radiances  already  compiled.  On  the 
other  hand,  from  the  Infrared  radiances  it  is  possible  to  draw  some  gen- 
eral conclusions.  In  the  infrared  it  is  not  unreasonable  to  treat  the 
lower  boundary  as  a blackbody  (emissivity  equals  unity).  Also,  for  a 
specified  altitude  "effective"  amounts  of  blackbody  clouds  can  be  intro- 
duced. The  effective  amount  is  actually  a product  of  the  amount  and  the 
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Beyond  the  scope  of  this  study 


emissivity,  although  apparent  emissivity  often  includes  subresolution 
cloud  holes.  Thus,  a 0.5  cover  of  cloud  with  a 0.5  emissivity  could  be 
associated  with  the  same  radiance  as  0.25  cover  of  blackbody  cloud  at  the 
same  altitude. 

All  of  the  radiances  listed  in  Table  6 can  be  manipulated  with  the 
above  assumption.  Tabulations  were  made  with  a blackbody  lower  boundary 
at  0,  3,  6,  and  9 km.  The  elevated  surfaces  might  be  used  to  represent 
low,  middle,  and  high  clouds.  Any  effective  amount  can  be  assigned.  For 
example,  if  0.5  cover  of  cloud  at  6 km  is  assumed,  the  "measured"  radiance 
would  be  equal  to  0.5  of  the  tabulated  results  for  6 km  and  0.5  of  the 
results  for  0 km.  The  sum  of  the  weighted  radiances  can  then  be  converted 
to  temperature  with  the  aid  of  Eq.  (25). 

When  this  type  of  analysis  is  pursued,  it  becomes  apparent  that  the 
infrared  channels  are  responsive  to  the  effects  of  high  clouds  but  un- 
responsive to  small  changes  near  the  lower  boundary.  Multispectral  tech- 
niques17* 30  that  have  been  applied  or  suggested  for  the  derivation  of 
cloud  height  or  amount  have  not  required  matched  spectral  responses.  The 
real  merit  of  matched  responses  is  in  application  to  the  detection  of 
known  spectral  variations  of  homogeneous  cloud  radiative  properties  es- 
tablished through  a proper  model.  The  only  other  advantage  is  that 
matched  sensors  will  respond  differently  to  a nonhomogeneous  scene  because 
of  the  spectral  variation  of  the  sensitivity  of  the  Planck  function  with 
temperature. 

As  an  example  of  the  use  that  can  be  made  of  data  in  the  radiance 
tables,  consider  the  spectral  regions  centered  at  455,  540,  and  1935  cm  ^ . 
The  regions  at  540  and  1935  cm  * have  already  been  examined  with  respect 
to  their  weighting  functions  in  winter.  From  the  midlatitude  winter  data 
new  radiances  (and  equivalent  temperatures)  were  formulated  from  the  tab- 
ulations to  correspond  to  0,  0.25,  0.5,  0.75,  and  1.0  equivalent  black- 
body cloud  cover  at  the  3,  6,  and  9 km  altitudes.  Figure  20  is  a plot  of 
the  Joint  radiances  (listed  in  Table  8)  for  the  455  cm  * and  540  cm  * 
channels.  The  heavy  dashed  curve  connects  joint  radiances  for  uniform 
complete  blackbody  coverage  at  each  of  the  indicated  altitudes.  Also 
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FIGURE  20  APPLICATION  OF  RADIANCES  FROM  TWO  INFRARED  CHANNELS  TO 
IDEALIZED  CLOUD  ANALYSIS  (MIDLATITUDE  SUMMER  ATMOSPHERE) 


» * 
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Table  8 
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SELECTED  RADIANCES  FOR  INDICATED  CLOUD  AMOUNTS 
AT  THREE  ALTITUDES 


-2  -1  -1  -1 

Radiances  (mW  m sr  (cm  ) ) 


Cloud 

Amount 

475  cm  1 

540  cm  * 

1935  cm" 

1 

3 km 

6 km 

9 km 

3 km 

6 km 

9 km 

3 km 

6 km 

9 km 

.00 

93.209 

93.209 

93.209 

103.188 

103.188 

103.188 

2.367 

2.367 

2.367 

.25 

92.704 

90.171 

86.222 

101.868 

97.439 

92.876 

2.235 

2.004 

1.869 

.50 

92.198 

87.133 

79.235 

100.547 

91.690 

82.563 

2.103 

1.643 

1.373 

.75 

91.692 

84.093 

72.246 

99.227 

85.940 

72.251 

1.971 

1.280 

0.875 

1.00 

91.186 

81.055 

65.259 

97.906 

80.191 

61.938 

1.838 

0.917 

0.377 

plotted  by  small  boxes  and  joined  with  straight  lines  are  the  joint  radi- 
ances at  6 km  and  9 km  for  each  of  the  cloud  amounts.  When  the  cloud 
height  is  held  fixed,  all  of  the  radiances  fall  on  a straight  line  con- 
necting the  joint  clear  column  radiance  and  the  joint  radiance  for  com- 
ple  te  coverage  at  the  given  altitude  (both  points  on  dashed  curve).  The 
fraction  of  distance  along  the  connecting  line  also  indicates  the  amount 
of  effective  cloud  cover.  (If  the  plots  are  made  in  equivalent  tempera- 
tures, then  the  points  do  not  divide  the  joining  line  by  portions  equiva- 
lent to  amounts.)  It  is  obvious  that  it  would  be  difficult  to  get  infor- 
mation at  low  altitudes,  and  that  the  ease  of  information  extraction 
increases  with  altitude.  The  plots  have  been  made  on  the  assumption  of 
single  cloud  layers  without  serious  noise  in  measurements. 

Whether  or  not  the  information  from  1935  cm-1,  for  example,  which 
nearly  matches  540  cm”^  in  response,  is  really  helpful  or  necessary  in 
dealing  with  amount  or  height  of  equivalent  cloud  has  not  been  ascertained. 
However,  it  is  possible  to  establish  some  difference  in  temperature  re- 
sponse to  a partial  cloud  cover. 

Since  the  Planck  function  is  much  more  sensitive  to  temperature  dif- 
ference at  higher  wave  numbers  (higher  frequencies)  than  at  low  wave  num- 
bers, the  warmer  (clear)  portions  of  a partly  cloudy  scene  will  contribute 
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relatively  more  to  the  total  radiance  at  1935  cm  than  at  540  cm  . In 

other  words,  the  Planck  function  for  the  1935  cm  interval  should  trans- 

late to  a higher  equivalent  temperature  over  a nonhomogeneous  scene  than 
that  for  the  540  cm  * interval,  even  with  a perfect  match  in  the  responses  * 

for  both  intervals  when  viewing  homogeneous  backgrounds.  Although  the 
responses  for  the  two  intervals  are  not  perfectly  matched,  the  concept  of 
the  differing  equivalent  temperature  can  be  demonstrated  from  the  computed 
radiances  in  Table  8.  If  TQ  represents  the  equivalent  temperature  for 
clear  skies  with  the  surface  at  0 km,  Tc  indicates  the  equivalent  temper- 
ature for  a complete  opaque  overcast  at  the  indicated  altitude,  and  or 
Tg  represents  the  equivalent  temperature  for  a partial  cloud  cover  at 
6 km  or  9 km,  then  the  following  temperature  ratios  are  obtained  for 
three  partial  cloud  covers: 


540  cm 


-l 


1935  cm 


-1 


(6  km)  (9  km) 


(6  km)  (9  km) 


Cloud 

Amount 


rvvi  r vii  rvvi  rvM 

lVTcJ  [VTcJ 


0.25 

0.50 

0.75 


.23 

.23 

.19 

.15 

.48 

.47 

.40 

.33 

.74 

.73 

.67 

.58 

262.9 

262.9 

T0  " 

265.1 

265.1 

243.3 

225.7 

T = 

243.1 

225.7 

A study  of  these  results  reveals  relatively  higher  equivalent  temper- 
atures at  1935  cm  * (in  comparison  to  540  cm  ) for  partial  cloud  covers, 


with  the  greatest  differences  for  the  highest  clouds.  Since  responses 


were  not  perfectly  matched,  the  equivalent  temperature  (Tq)  for  1935  cm 
over  a homogenous  surface  under  clear  skies  was  2.2k  higher  than  for  540 


-1 


cm  . With  an  opaque  uniform  cloud  cover  at  9 km  the  equivalent  temper- 


atures (T  ) for  the  overcast  state  were  the  same.  When  partial  cloud 
c _1 

covers  are  introduced  at  9 km,  it  is  apparent  that  at  1935  cm  only  157. 
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of  the  total  clear-overcast  temperature  decrease  occurs  by  257.  of  the 
possible  cloud  cover,  and  only  337.  of  the  temperature  decrease  occurs  as 
the  cloud  cover  Is  Increased  to  50  percent.  Comparable  figures  for  540  cm 
are  237.  and  477.  of  the  total  temperature  decrease.  Since  the  percentage 
temperature  decline  is  smaller  at  1935  cm  * for  each  quartile  of  cloud 
amount,  the  temperatures  are  remaining  relatively  higher.  The  excess  of 
equivalent  temperatur*  at  1935  cm  ^ over  that  at  540  cm  ^ increases  from 
2.2K  for  clear  skies  to  a maximum  of  6.7K  for  50  percent  cloud  cover. 
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VIII  CONCLUSIONS 

This  study  has  demonstrated  the  feasibility  of  modeling  simplified 
but  reliable  transmittance  representations,  over  infrared  spectral  inter- 
vals on  the  order  of  20  wavenumbers,  for  applications  to  the  nonhomogen- 
eous  molecular  atmosphere.  Parameterizations,  which  include  temperature- 
dependent  scaling  coefficients,  are  based  on  detailed  monochromatic 
absorption  coefficients  computed  from  pertinent  absorption-line  parameters 
at  a variety  of  pressures  and  temperatures.  Data  and  the  transmittance- 
radiance  models  are  available  for  extension,  refinement,  and  empirical 
testing  as  well  as  for  introduction  into  a meaningful  software  system 
package  for  future  design  studies. 


The  transmittance  representations  were  used  in  the  generation  of 
background  information  that  documents  radiance  variations  with  changes  in 
spectral  position,  atmospheric  structure,  viewing  angle,  and  boundary 
conditions.  Infrared  radiances  and  their  equivalent  temperatures  that 
were  simulated  in  the  presence  of  opaque  surfaces  at  various  altitudes 
may  be  examined  to  infer  gross  effects  of  cloud  height  and  amount.  The 
higher  the  effective  opaque  surface,  the  more  responsive  were  the  infrared 
sensors.  The  microwave  brightness  temperatures  and  their  limb  variations 
over  surfaces  with  reduced  emissivity  were  dependent  on  the  type  of  re- 
flecting surface.  Limb  brightening  was  reduced  (darkening  enhanced)  over 
diffuse  surfaces  relative  to  specular  surfaces. 


Considerable  emphasis  was  given  the  identification  of  spectral  inter- 
vals at  widely  separated  wavenumbers  with  the  same  general  response  to  the 
cloud-free  atmosphere.  Although  some  very  good  matches  were  obtained,  it 
is  difficult  to  achieve  a meaningful  match  unless  both  regions  have  the 
same  distributions  of  absorber  types.  Without  such  similarity,  a match 
found  under  one  condition  would  not  hold  for  another.  For  those  interval 
pairs  with  responses  from  near  the  lower  boundary,  the  spectral  variation 
of  the  continuum  attenuation  causes  a matching  problem.  When  a match  is 


found,  the  spectral  variation  of  the  temperature  sensitivity  in  the 
Planck  function  causes  the  higher  wavenumbers  to  show  higher  temperatures 
over  a nonhomogeneous  scene.  Different  intervals  with  matching  cloud- 
free  responses  will  be  useful  for  the  identification  of  spectrally  vary- 
ing radiative  properties  of  clouds,  once  proper  physical  cloud  models 
are  introduced.  High  priority  should  be  given  to  the  application  of  a 
suitable  cloud  model. 
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Appendix  A 


Program 

RDGOAT 

DAVISV 

COMB IN 

VTRAN 


TAUFIT 


RADIANT 


SUMMARY  OF  COMPUTER  PROGRAMS 


Function 


Reads  AFGL  tape  of  absorption  line  parameters,  ex- 
tracts parameters  for  materials  and  wavenumber  of 
interest.  Computes  ratios  of  Lorentz  to  Doppler 
lines,  lists  and  stores  data. 

Generates  extensive  files  of  monochromatic  absorp- 
tion coefficients  for  specified  absorbers,  tempera- 
tures, pressures  and  spectral  spacing.  Stored  on 
tapes . 

Weights  absorption  coefficients  for  uniformly  mixed 
gases  with  ratio  of  their  mixing  ratio  to  that  for 
CC^.  Sums  ..sighted  coefficients  within  same  spectral 
region  and  restores  in  single  file. 

Uses  absorption  coefficients  from  DAVISV  or  COMBIN  in 
conjunction  with  specified  absorber  amounts  to  pro- 
duce spectrally  averaged  transmittance  (monochromatic 
transmittances  convolved  with  any  stipulated  spectral 
response  function).  Lists  and  cardpunches  final  trans- 
mittances. 

Two  routines  that  apply  a nonlinear  least  squares  curve- 
fitting routine  to  data  from  VTRAN.  TAUFITl  applied  at 
single  reference  temperature;  solves  for  all  parameters. 
TAUFIT2  holds  one  parameter  fixed  and  solves  for  the 
other  parameters  at  all  temperatures,  and  computes  tem- 
perature scaling  coefficients  to  complete  model. 

In  application  to  arbitrary  atmospheric  structures,  uses 
modeled  transmittance  representation  with  coefficients 
for  recovery  of  transmittance  from  transmittance  func- 
tion. Parameters  for  computing  continuum  transmittance 
also  are  provided  along  with  any  other  special  infor- 
mation, such  as  slant  path  desired.  Transmittance 
(including  continuum)  and  weighting  function  profiles 
generated;  radiances  computed. 
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SUMMARY  OF  COMPUTER  PROGRAMS  (concluded) 

Program  function 

ABS02  0 (5  mm)  microwave  absorption  coefficients,  trans- 

mlttances,  weighting  functions,  brightness  temper- 
atures . 

ABH20  Same  function  but  for  H.^0  in  vicinity  of  either  22 

GHz  or  183  GHz.  Modified  version  (GAMAC)  provides  a 
card  deck  of  absorption  coefficients  for  combination 
with  0^  coefficients. 

M1CR0W  Combines  0^  and  H (0  absorption  coefficients  to  provide 
joint  transmittances. 

PLOTAB  Plots  absorption  coefficients. 

TAUPLT  Plots  data  points  and  corresponding  transmittance  repre- 
sentation. 
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RECOVERY  OF  TRANSMITTANCE  FROM  THE 
TRANSMITTANCE  FUNCTION 


The  transmittance  function  Y is  defined  by  Eq.  (12)  in  the  text. 

c 

To  avoid  table  look-up  and  interpolation,  a numerical  routine,  based  on 
a Chebyshev  polynomial  expansion  available  in  FORTRAN  (IBM  Scientific 
Subroutine  Package),  is  introduced  for  rapid  recovery.  Seven  coefficients 
of  the  Chebyshev  polynomial  expansion  appeared  to  be  sufficient  to  ensure 
at  least  three-place  accuracy  in  t for  the  worst  case.  A set  of  coeffi- 
cients, along  with  a pair  of  scaling  coefficients  xQ  and  x^  are  required 
for  each  value  of  parameter  C.  In  practice  the  logarithm  of  the  trans- 
mittance function  is  determined,  so  that  the  coefficients  were  computed 

to  relate  £n(-ini)  to  £nY  . 

c 

The  transmittance  function,  specified  either  by  Eq.  (13)  or  Eq.  (14), 
is  determined  from  atmospheric  data  and  a known  parameter  C.  The  recovery 
procedure  begins  by  defining 

E = 2(x.£nY  + x.) 

1 c 0 

and  setting  H^H^O.O,  initially.  Then  the  quantities  H2  and  HQ  are  de- 
termined from  the  following  computational  loop  from  i=l  through  i=N: 

H2  " H1  i H1  = H0 

j = N + 1 - i 
H0  - EHj_  - «2  + C(j) 


•h<  i C(  ) Is  the  set  of  Chebyshev  coefficients  (j  = l,  2 N).  After 

- nit  run,  the  resultant  »2  and  HQ  determine  the  quantity  F 
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F s An(-£nT)  = 0.5^C(1)  - H2  + HQ  j 

and  T = exp ( -exp  F) . 

Two  special  cases  arise.  For  the  parameter  C=0.0 

r = exp(-Y0) 

and  for  C=1.0 

T = 1 - erf  C/5T  Y/2) 

Table  B-l  lists  all  of  the  coefficients  of  the  Chebyshev  polynomial 
expansions  and  required  scaling  coefficients.  These  coefficients  relate 
only  to  the  transmittance  function,  not  the  atmospheric  representation. 
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COEFFICIENTS  FOR  NUMERICAL  RECOVERY  OF  TRANSMITTANCE  FROM 
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APPENDIX  C 

ADDITIONAL  WEIGHTING  FUNCTION  PROFILES 
COMPUTED  FOR  SELECTED  OXYGEN  AND 
WATER  VAPOR  MICROWAVE  FREQUENCIES 


HEIGHT  — km  HEIGHT 


